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Some ten years ago only the few were aware of the 
existence of molecules without chemical bonds. 
Nowadays a new scientific discipline, chemical to- 
pology, appeared to study them. This science has suc- 
ceeded in predicting and producing many remarkable 
chemical compounds (e.g. catenanes whose molecules 
are interconnected like links of a chain), in denning 
the electronic structure of molecules more precisely, 
and so on. 

The complicated concepts of chemical topology 
are presented here in a clear and comprehensible 
form. ' 
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sity and high-school students majoring in chemistry 
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"The more I study . . . the more I become 
impressed with the harmony ... which 
exists between the chemical and al- 
gebraical theories. ... There is an 
untold treasure of hoarded algebraical 
wealth potentially contained in the 
results achieved by the patient and 
long-continued labour of our uncon- 
scious and unsuspected chemical fel- 
low workers." 

J. J. Sylvester 
Amer. J. Math., Vol. 1, 64 (1878) 



Author's preface 



In recent chemical literature more and more often one 
comes across such terms as "topology of a molecule", 
"topological properties", "topological bonding", etc. Broad- 
ly speaking, topology is a branch of mathematics deal- 
ing with the phenomenon of continuity. A more detailed 
definition would require from the reader a more pro- 
found knowledge of many complex mathematical concepts. 
This book, however, is designed primarily for a chemist. 
Besides, it is a popular science book (of course, if one 
regards as popular science not only books of what-chemis- 
try-makes-of -timber type). For this reason we shall give 
here only the most essential data on topology. 

High-school geometry, as a rule, looks into those prop- 
erties of geometrical figures that are associated with 
the notions of length, angle, area, and volume. These 
properties of figures are called metric. In addition, there 
are properties of figures that are independent of their size 
or shape; for instance, stretching, shrinking,$bending, 
and twisting neither bring various points of the figure 
into contact, nor induce any ruptures in it. One calls 
such transformations homeomorphic. The transformation 
of a doughnut into a mug shown in Fig. 1 (p. 8) is homeo- 
morphic. A juggler in Fig. 2 (p. 9) removing the vest 
without taking off the jacket also performs a homeomor- 



phic transformation with the vest. The properties of 
bodies remaining invariable in homeomorphic transfor- 
mations are examined by topology and called topological. 
The shape and the size of a figure are of no importance 
in topology. 

Metric parameters of molecules, such as equilibrium 
distances between nuclei, valence angles, etc., are closely 
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connected with their electronic structure. As will be shown 
later, however, some molecular properties depend weak- 
ly, if at all, on geometry and are determined essentially 
by the way the atoms are connected to one another, i.e. by 
topology. For example, both hydrogen atoms in a water 
molecule form covalent bonds with the oxygen atom and 
are not bonded together. 

An important point should be raised here, though. 
Strictly speaking, the term "topology of a molecule", as 
well as other related terms, cannot be regarded as wholly 
appropriate since chemical structures are discrete, and it 
would be more competent to consider a molecule in terms 
of the graph theory. Surely, a graph can be defined as 
a topological object, but the methods of the graph theory 
are not topological for the most part. Nevertheless, chem- 
ists often use the above-mentioned, not very rigorous 
terms not only out of the vainglorious urge to show off, 
but primarily to stress the fact that some properties of 
a molecule are independent of its geometry and the nature 
of the chemical bond. In the same sense, we shall also be 
dealing here with chemical topology. 

It should be noted that this book, both in its structure 
and in its content, is in many respects logically related 
to our earlier booklet: I. S. Dmitriev, Symmetry in the 
World of Molecules, Mir Publishers, Moscow, 1979. This 
relation is not accidental since both topology and symmet- 
ry of molecules furnish a valuable qualitative and semi- 
quantitative complementary information about their 
structure and properties. 



8 




Fig. 2 



Chapter One Approximation f ' 

to Accuracy 
and Accuracy 
of Approximation 

The Value of a Decimal Place 

Modern quantum chemistry has developed diverse 
methods of calculating molecules. Usually these meth- 
ods are divided into three groups: 

(1) comprehensive nonempirical ab initio calculations 
involving all electrons of a system; only the values of 
internuclear distances, if at all, are borrowed from ex- 
perimental data; 

(2) semi-empirical semiquantitative calculations in 
which only valent electrons (and sometimes even not all 
of them) are taken into account; in this case some of the 
integrals are discarded and others are approximated by 
a special selection of parameters whose values are bor- 
rowed from experimental data; 

(3) qualitative methods of analysis based on^ the fun- 
damental concepts of quantum chemistry, the symmetry 
considerations, and, more and more often recently, on 
the topological methods. 

Surely, all three approaches to the problem of elec- 
tronic structureTof molecules?are closely interrelated. The 
modern theory of chemical bonding, on the one hand, 
relies on the growing utilization of powerful computers 
enabling various systems offinterest to a chemist^to be 
calculated in greater detail, and, on the other "hand, 
develops newjqualitative methods, thus enriching itself 
ideologically. 

We begin with a brief outline of some methods and 
approximations of quantum chemistry. 

The LCAO MO method. The method of molecular orbitals 
(MO) dominates in modern quantum chemistry. The main 
idea of the method is that the electrons in a molecule 
are accommodated in definite MOs just as in an atom they 
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are accommodated in definite AOs, i.e. in accordance^, 
with the Pauli principle, the orbitals being usually 
filled as the energy of the MO grows. 

Normally MOs are^made up as linear combinations of 
atomic orbitals (LCAO) of atoms composing the system, 
i.e. are written in the form 

*|>i = 23 c,iq>j, (1) 
j=i 

where i is the number of the MO i|v, / are the numbers of 
atomic q>-orbitals; c t j are the numerical coefficients defin- 
ing the contributions of individual AOs into the giv- 
en'MO. 

Such ajway ©^constructing a MO is based on the assump- 
tion that an atom represented as a definite set of orbitals 
remains distinctive in the molecule. 

The complete iV-electron wave function of a molecule 
is represented in the form of the Slater determinant com- 
posed of the MOs: 

Y (x u . . . , x N ) — 

= \/\/W\ 



^1 Ifl (*2> 

1|>2 fa) $2 



^1 (*n) 

^2 (X N ) 



\p N (Xi) \lf N (x 2 ) ... \p N (x N ) 



(2) 



Electrons interchanged. When the many-electron func- 
tion of a molecule is written in the form of a determinant, 
the fundamental antisymmetry principle (the Pauli exclu- 
sion principle) of quantum mechanics is satisfied. Accor- 
ding to that principle an A^-electron function must be 
antisymmetric, i.e. it must change sign whenever spatial 
and spin variables of any two electrons are interchanged: 

^¥ (xi, x 2 , . . ., Xi, . . ., Xk, . . ., Xff) = 

= W (xi, x%, . • ., £jt, • • ., Xi, . . ., Xff). 

This is exactly the property that expression (2) possesses 
since any determinant changes its sign when its two 
rows, or two columns, are interchanged. (The readers 
themselves may easily check this in the case of a deter- 
minant of the second order.) It thus follows that each 
molecular spin orbital can accommodate not more than 
one electron and accordingly each molecular orbital can 
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accommodate not more than two electrons (witt* opposite 
spins). Otherwise, two rows in determinant (2) become 
identical, and such a determinant with two identical rows 
(or columns) is always equal to zero. And what if the 
Slater determinant is equal to zero? This means that 
*F x 2 , ■ ■ ., x N )= 0, i.e. the system in question can- 
not exist in such a state. 

Thus we came to the conclusion that in the MO method 
the total (JV-electron) wave function of a molecule is 
expressed via one-electron functions, that is, molecular 
spin orbitals which, in their turn, are expressed via atom- 
ic spin orbitals. Subsequently we shall be utilizing only 
MOs expressed in the LCAO form. Although the spin 
variables will not be taken into account, the Pauli prin- 
ciple will be obeyed. 

In 1928 the British scientist D. Hartree suggested his 
method of calculating orbitals. In 1929 his approach was 
generalized by V. A. Fock who derived more general equa- 
tions, putting to use the principle of antisymmetry. 

Roothaan's equations. In 1951 using Fock's method the 
American physicist C. C. Roothaan worked out a system 
of nonlinear algebraic equations providing the AO coeffi- 
cients of Eq. (1): 

2(^-8^)^ = 0, (3) 



where F ab is the quantity specifying the contributions 
that various interactions (electron-electron, electron- 
nucleus, etc.) make to the energy of a molecule; since 
this quantity has a fairly complicated mathematical form, 
it is not presented here in an explicit form; S ab is the 
orbital overlap integral; e ft is the energy of the k-th MO. 
Equations of type (3) are called secular. 

Let us compose the following column matrix from the 
coefficients of expansion of the i-th MO into the corre- 
sponding AOs: 



/'a \ 



Then Sjp. (1) can be written in a compact form as 

where <p is the row of the AO: 

<P = (<Pi» <Pa> • • m <Pm)- 

Assembling the columns c t into the matrix C of m X n 
size (m is the number of AOs and n is the number of 
MOs), we can rewrite Eq. (1) in the form of a single mat- 
rix equality 

yjp -= <pC, 
where i|> is the row of the MO: 

Then the Roothaan equations can be written in the 
matrix form: 

FC - SCe (4) 

where F is the so-called Fock matrix composed of the 
elements F ab ; S is the overlap matrix; e is the matrix 
whose diagonal elements specify the MO energies, all other 
elements being equal to zero. 

Solution of the Roothaan equations calls for laborious 
computations which become more and more so as the ba- 
sis broadens, i.e. the number of AOs increases. In the 
last two decades theoreticians have invented a great num- 
ber of clever tricks in their attempts to find how to cal- 
culate molecules "the Roothaan way". However, the 
system of equations (4) turned out to be very tough to 
handle. Despite the ingenuity of researchers and the 
advancement in computer technique practical application 
of the Roothaan method is substantially limited by the 
size of molecular systems. 

Finding a way out. Calculations can be considerably 
simplified by using approximation of the zero differential 
overlap (ZDO). The essence of this approximation is that 
the overlap of different AOs i|5j and tyj is assumed to be 
equal to zero for any element dx of the volume of a mole- 
cule: 

q>i<Pjd% = (i 
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Consequently, all overlap integrals also tur/ to zeroi 

^ ii = \ WtWjdT; = 

as well as many other integrals, including those which are 
the hardest to crack. Initially, the ZDO approximation 
was introduced as a postulate, and only later was it 
substantiated theoretically and its validity limits were 
set. Many varieties of the ZDO method are now extensively 
employed in quantum chemistry. 

Cost of accuracy. At present with the appearance of the 
ZDO approximation that considerably simplifies calcula- 
tions, theoreticians often prefer a semi-empirical approach 
in which most of the "troublesome" integrals are either 
neglected completely or expressed via the parameters found 
experimentally (orbital ionization potentials, elec- 
tron affinity etc.). In this case the successful calibration 
of the empirical parameters can offset the loss in accuracy 
caused by various simplifications introduced into the 
Roothaan calculation method. 

Today semi-empirical calculation methods provide 
a great deal of theoretical information about the electron- 
ic structure of chemical compounds. This preference for 
semi-empirical methods is primarily due to the following 
two factors. 

First, nonempirical calculations consume much compu- 
ter time which makes them more expensive. 

Second, semi-empirical methods admit, as a rule, more 
readily, descriptive interpretation in terms of conven- 
tional chemistry. 

While giving credit to these methods one should not 
forget their drawbacks, in particular, the lack of rigid 
criteria to assess their accuracy. According to the British 
scientist G. G. Hall, each theoretician has a semi-empiri- 
cal method he believes in, but never trusts other semi- 
empirical methods. 

Complexities of a Simple Method 

The first and the simplest form of the molecular orbit- 
al theory was proposed by E. Hiickel in 1931. Despite its 
simplicity (which is largely seeming, anyway) this meth- 
od is still extensively used in theoretical chemistry to- 
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day in it^ various modifications, emphasizing once again 
the fact that the calculation methods cannot be uncondi- 
tionally judged as the good ones or the poor ones for 
much depends on what is the system, which of its prop- 
erties are studied, and what approximation is used. 
In quantum chemistry, as in other disciplines for that 
matter, it is important not only to attain the maximum 
possible accuracy of calculations, but also to choose the 
approximation most adequate for the given problem. 
It is relevant here to recall the words of the eminent 
French mathematician A. Poincare: "The problem is 
not in what is the answer, the problem is in what is the 
question." 

Because of its simplicity and limited computational 
requirements the Hiickel theory was especially useful 
in the period before the application of computers to quan- 
tum-chemical problems (1930-midl950s). However, 
since the 1960s the two opposite views concerning the 
Hiickel theory have been advocated. One is that the days 
of this theory are numbered because high-speed computers 
and more sophisticated MO methods are now available. 
The other view is that there is still room for using the 
Hiickel theory in organic chemistry, especially on a qual- 
itative level as a guide for the chemists in planning and 
interpreting experiments. This latter opinion is consis- 
tent with the view that most experimental chemists pre- 
fer a "pencil-and-paper" method for their everyday re- 
search. 

The Hiickel method. In the Hiickel molecular orbital 
method the Fock matrix (see Eq. (4)) is replaced by 
a certain "effective" matrix H whose elements H hj are not 
explicitly specified. 

Then each MO is composed of TV valence AOs of atoms 
in a molecule, i.e. is written in the LCAO MO form. 

Frequently the secular equation defining the coefficients 
c i} and the orbital energies e t is also derived with the use 
of the ZDO approximation. Then the equations of the 
Hiickel method take the following form: 

sVw-8«6 w )c,, = (5) 

3=1 

(4=1,2, ...,N), 



where 



0, if k^j, 



f 0, if k^j 
6w= l 1, if * = /. 



Eqs. (5) have nontrivial solutions provided the secular 
determinant is equal to zero: 

det | H h j — Eib hJ | = 

(ff«-e,) 



= det 



H 



Zi 



H i2 

(# 22 — 8j) ... 



H 



2N 



H 



Nl 



= 0. 



(6) 



H N2 ■•• (Ann — e s) 

As has been mentioned, the matrix elements H h} are 
not specified explicitly in the HMO method. They are 
treated as parameters whose values are determined from 
experimental data. 

The same equations in a matrix form. Eqs. (5) and (6) 
can be written in a more compact form if one uses matrix 
notation. Thus, Eq. (5) in matrix form is 

(H - e| I)c, =0, (7) 

where Cj is a column matrix composed of the AO coeffi- 
cients; H is a matrix composed of the elements H k y, I is 
the diagonal unitary matrix. 

Eq. (6) takes the following form in matrix notation: 



det | H - z t l | = 0. 



(8) 



Overlap is allowed for. In the Hiickel method the overlap 
integrals are usually discarded in calculations. But they 
can be taken into account. In this case Eqs. (7) and (8) 
take a somewhat complicated form: 

(H - ei S)c e = 0, (9) 

det | H - e;S | = 0. 

However, account taken of the overlap integrals does 
not generally improve the calculations appreciably. 

Maximum simplicity. So let us return to the simplest 
version of the Huckel method. We have not as yet men- 
tioned still another simplifying assumption of this meth- 
od: most of the matrix elements in Eq. (6) are assumed 
""*' "-" *«.. 
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to be eq^l to zero; the remaining elements are the diag^. 
onal elements H n called the Coulomb integrals and deuo- 
ted by a; and those nondiagonal elements H l} (i =fcj) 
for which the indices i and / relate to the orbitals of the 
chemically bonded nearest-neighbour atoms. The nondiago- 
nal matrix elements H i} are called the resonance integrals 
and denoted by the symbol p^. 

The Coulomb integral a t defines the AO <pi and, con- 
sequently, the atom to which this AO belongs; in this 
case a < 0. The resonance integral defines the AOs <p { 
and <py, in the case of finite internuclear distances R 

In the HMO method the total energy of a molecule in 
the ground state is equal to the double sum of orbital 
energies: 

£^ = 2^, (10) 
n 

(over occupied MOs) 

The factor 2 reflects the fact that we consider a molec- 
ular system with closed electronic shell, that is, each 
MO accommodates two electrons with opposite spins. 

Eq. (10) is obviously erroneous because in the one- 
electron approximation the total energy does not reduce, 
in fact, to the sum of the orbital energies, i.e. 



but differs from that sum by the value of the averaged in- 
terelectronic repulsion energy. The indicated drawback 
of the HMO method, as well as its other weak points, 
naturally limits the application of this method for the 
quantitative analysis of electronic structure of molecules, 
but it is quite acceptable for qualitative conclusions. 

The HMO method proved to be most efficient in studies 
of conjugated and aromatic systems. 

Properties of conjugated systems. Conjugation manifests 
itself in molecules in many ways. By formal definition, 
conjugated compounds are those which feature alternating 
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singleand multiple bonds, e. g. butadiene/ CH 2 = 



=CH— CH=CH., or benzene 



As a rule, mole- 



cules of this kind exhibit a number of peculiar properties, 
primarily, higher thermodynamic stability. How much 
"higher"? Comparing to what? When higher stability of 
conjugated systems is discussed, chemists compare them 
to hypothetical systems possessing noninterfering localiz- 
ed multiple bonds. For instance, the heat of hydration of 
one double bond in a "real" butadiene molecule is less 
than the heat of hydration of ethylene, and "the heat of 
combustion of butadiene is less than that of the hypo- 
thetical system with localized double bonds by apprecia- 
ble (although small) value of 14.7 kJ/mol. In the case of 
benzene the respective difference is much higher and 
amounts to 155 kJ/mol. 

Thus, the first characteristic of conjugated system is 
associated with their energy parameters, e.g. the energy 
of butadiene and benzene is lower than that of the re- 
spective classical structures. 

The second characteristic of conjugated molecules is 
their inherent tendency to equalize the lengths of sin- 
gle and multiple bonds. For instance, whereas in 

ethane R c _ c = 1.54 A and in ethylene R c =c = l-34 A, 
the interatomic distances in butadiene are as shown 
below: 

1.35 A 1.48 A 1.35 A 

H2C=^=^CH CH CH2 

In benzene the lengths of all carbon-carbon bonds are 
equal to 1.397 A. 

Other characteristics of conjugated molecules can also 
be indicated: the shift of the UV absorption maximum of 
butadiene and other polyenes to the longwave region of 
spectrum as compared to ethylene, the lowering of ioniza- 
tion potentials, and so on. In short, .there is much evid- 
ence pointing to the fact that the properties of conjugated 
systems are determined by the interaction of multiple 
bonds. 

However, as will be shown in Chapter 4, not all formal- 
ly conjugated systems possess the above-mentioned prop- 
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erties, stN, that the definition of the term "conjugation^ 
becomes rather difficult. 

Note that most of theoretical concepts in chemistry are 
associated with a certain approximation and each approx- 
imation, in its turn, with the introduction of new notions 
whose physical meaning is often far from being clear. 
Among such notions are electronegativity, aromaticity, 
conjugation, hybridization, and others. Experimental 
chemists impart significant empirical meaning to those 
notions, and theoretical chemists manage to interpret 
various aspects of that meaning, more or less convincing- 
ly, using a certain model. 

The Hiickel approximation. One of the most popular 
models is the model of unsaturated molecules which is 
based on the Hiickel approximation*, or the approximation 
of CT-n-separability. The presence in conjugated systems 
of a symmetry plane coinciding with the plane of the 
nuclear skeleton of a molecule allows its MOs to be divided 
into two orthogonal groups, symmetrical (a-MOs) and 
antisymmetrical (ji-MOs) relative to the reflection in that 
plane. Consequently, one can consider the n shell of 
a molecule separately from the a frame. In other words, 
in the jt-electron approximation we consider not just 
valence orbitals but valence n orbitals. 

In this case the a frame is assumed to be conservative, 
in the sense that it is not polarized and its orbital energies 
are lower than those of the ji shell. However, the non- 
empirical calculations demonstrated that the a frame is 
not so much conservative as our ideas about that frame 
are. The a electrons contribute substantially to the 
energy and appreciably to the polarizability of conjugated 
systems. Nevertheless, the jt-electron approximation is 
quite acceptable in those numerous cases when qualita- 
tive and semi-quantitative results must be obtained. We 
should only remember that what we use is an approxima- 
tion! 



* The Hiickel approximation (not to be confused with the 
computational Hiickel method I) is also called the ji-electron appro- 
ximation. 
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Calculations Without a Computer / 

Let us consider now a few examples of calculations of 
conjugated hydrocarbons by the Hiickel method and in 
the Hiickel approximation. 

The first example: ethylene. In an ethylene molecule 
there are two ji electrons forming one n bond between the 
carbon atoms. The corresponding ji-MO can be represented 
as a linear combination of two 2p n AOs of carbon atoms, 
q>! and cp 2 : 

■ti = Cii<p t + c 2i <p 2 - 

To write the secular equation, one has to label the 
carbon atoms in an ethylene molecule: 

I 2 
C-C 

and then to compose the determinant whose diagonal ele- 
ments are all equal to (a — e) and nondiagonal elements 
located at the intersection of the i-th row and the j-th 
column are equal to either p or depending on whether the 
i-th and /-th atoms are chemically bonded (via a single 
or a multiple bond) or not. Then the following secular 
equation can be ascribed to ethylene: 

a — e p 

P a — e 

from which, by solving the determinant, we obtain: 
( a _ e )« _ p2 = 0. 

This equation has two roots: 
e x = a + p, 
e 2 = a — p. 

Now we can find the coefficients c lt and c 2i . In our case 
Eq. (7) takes the form 

c«(a — e i ) + c 2i P = 0. 

Substituting into this equation first e x and then e 2 
we get 

for ej Cu = c 12 = c and \p t = c (ip, + q> 2 ), 
for e 2 c 21 = — c 22 and i|) 2 = c (q>( — (p 2 ) . 



det 



= 0, 



20 



Then taking account of the condition for normalization 

of orbitals, we obtain c — l/J/2. 

Thus, we finally get the following expressions for 
n-orbital energies and for n MOs of ethylene: 

8 1 = a + p Vl = 1/1/2(9! + q> 2 ), 

and 

e 2 = a — P i|j 2 = 1/1/2 (q)! — q> 2 ). 

Since both parameters a and {$ are negative, the molec- 
ular orbital ^ occupies a lower energy state than the 
orbital i|) 2 : e x < e 2 . 

Figure 3 illustratesjboth MOs together with their squares. 
It can be seen from this figure Jthat the MO ^ is char- 




Fig. 3 



acterized by a high electronic density between the nuclei 
whereas the MO t|) 2 , and its square, turn into zero at a point 
located between the nuclei. Accordingly, the n MO ifo 
is called the bonding (BMO) and the n MO a|) 2 antibonding 
orbital (ABMO). 

The second example: butadiene. In the case of a butadi- 
ene molecule 



1 2 
H.C=CH- 



3 4 
-CH=CH 2 



the HMO method provides the following secular equation 

a— e p 
P a— e P 

P a ~ e P 
p <* — e 



det 



= 0. 



(11) 



For the sake of convenience we shall introduce the 
designation 



a—e 



P 



Then Eq. (11) takes the form 
i 100 



det 



1x10 
1x1 
1 x 



:0. 



Solving this determinant, we easily obtain the biquad- 
ratic equation 

x 4 — 3x 2 + 1 - 

whose roots are 

t±/5 



^1,2 — 



_ 1 + /5 



Correspondingly, we obtain four energy levels (in the 
order of their increase): 



e, = a+1.618B 
e 2 = a + 0.618B 
e 3 = a — 0.6180 
e 4 = a— 1.6186 



| energies of bonding ji MOs, 
J energies of antibonding n MOs. 



Spending some more time, we can obtain expressions 
for the corresponding four MOs*: 

^ = 0.3717 (p x + 0.6015 (p 2 + 0.6015 <p 3 + 0.3717 q> 4 , 
^ 2 = 0.6015 + 0.3717 (p 2 -0.3717 <j> 3 -0.6015 <p 4 , 
t|) 3 = 0.6015 -0.3717 <p 2 -0,3717 q> 3 + 0.6015 q> 4 , 
% = 0.3717 q)j —0.6015 cp 2 + 0.6015 (p 3 — 0.3717 q> 4 . 



* If a molecule possesses certain elements of symmetry, MOs 
are easier to calculate and classify using the group theory (see 
our book Symmetry in the World of Molecules). 
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From Tig. 4 it can be seen that the it orbital hSe-no 
nodes and is strongly bonding, and the n MO tp 2 has 0*C&- 
node. The latter orbital bonds the atom C 1 with G 2 and " 
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the atom C 3 with C 4 loosening the bonding between the 
atoms C 2 and C 3 . The bonding effected by this MO is 
weaker compared to that effected by the former MO. The 
it MO \|) 3 is a weak antibonding orbital which strengthens 
only the C 2 — C 3 bond. And finally, the antibonding ji MO 
o|) 4 loosens all the three bonds, 



Derealization energy. To describe quantitatively the 
interaction of it bonds in conjugated systems, the concept 
of derealization energy (resonance) E n is frequently 
introduced. Usually it is denned as follows. 

The derealization energy E D is the difference 
between the total n-electron energy of a^given con- 
jugated compound (Efotai) and the energy of all 
localized bonds in that compound (Eu>m)*- 
For example, in butadiene four n electrons are accom- 
modated in the two bonding MOs, ^ and ip 2 , and the 
total n-electron energy of the molecule is 

Emm = 2 (a 4 1 .6180) + 2 (a + 0.6180) = 4a + 4.4720. 

Assuming the jt bonds to be noninteracting, we find 
that EmTi is equal to the double ji-electron energy of 
ethylene, i.e. 

£foS = 4a + 4p. 

Then the derealization energy equals 0.4720. Since 
< 0, the derealization results in the lowering of the 
total energy of a butadiene molecule, i.e. increases its 
stability. 

Polyene with N carbon atoms. In the general case of 
polyene with N carbon atoms forming a linear carbon 
chain [CH 2 (CH) JV _ 2 CH 2 1 the secular equation takes the 
form 

x 1 ... 
1 x 1 . .. 



... 1 
... x 
and its solution is given by the general formula 

(fi = l, 2, ...,N) 

* Somewhat different definitions that are sometimes given will 
not be discussed here. 
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or 



= a + 2pco8 



The AO coefficients in the expression for the k-th MO 
are found from the general formula 



z Hh : 



r 



sin 



\ N+i I ■ 



When the chain is closed. Using the Hiickel method and 
the it-electron approximation, let us consider now the 
most typical aromatic compound, benzene. Repeating the 
arguments and the calculations presented above, we ob- 
tain the secular equation for benzene in the form 



det 



= 0. 



Whence: 



H 1 
1x10 
1x10 
1x10 

1 x1 

1 1 x 



(x 2 -l) 2 (x 2 — 4) = 



and 



Xj, 2 — i 2, x 3 4, 5, e — i 1 • 
Correspondingly, we get six energy values of jt-MOs*: 

e,=-~a + 2P; e 2 = e 3 -a + p; 

e 4 = e 5 = a — p and e 6 = a- 



-2P. 



The six n electrons in a benzene molecule in the ground 
state occupy the levels e l5 e 2 , and e 3 , so that Emai 
(benzene) = 6oH 8p. This value differs by 2p from the 
n-electron energy (6a + 6P) of a system with three loca- 
lized jt bonds. Therefore, E D (benzene) = 2p. 



* The sequence of e's subindices corresponds to the growth of 
the values of e and does not coincide with the sequence of x's 
subindices. 



Quite often, when calculating the derealization ener- 
gies of cyclic polyene systems, the latter are compared with 
isoconjugated acyclic polyenes, e.g. benzene with hexa- 
triene etc. 

In the general case the following expressions for orbital 
energies of monocyclic polyenes can be obtained: 

£(1 = a + 2pcos(-^i) (12) 

(u. = l, 2, .... AT). 

Substituting the indicated values of u. into Eq. (12), 
we get the n-orbital energies of the system. 

Cyclobutadiene. For a cyclobutadiene molecule (N = 4) 
Eq. (12) yields the following results: the molecule con- 
tains one bonding n-MO (e x = a + 2EJ) accommodating 
two electrons, two nonbonding rc-MOs (e 2 = e 3 = a) 
accommodating one electron each, and one antibonding 
vacant n-MO (e 4 = a — 2|i). 

Thus, cyclobutadiene must be a very unstable and 
reactive compound, as it is indeed the case.* 

Willstatter's discovery. In 1905 the German chemist 
R. Willstatter obtained from the pseudopelletierine 
alkaloid a substance to which he ascribed the following 
structure: 




This is cyclooctatetraene. However, Willstatter was 
distrusted because it was considered, in accordance with 
the conceptions of his time, that a cyclic polyene must 
exhibit the "benzene-like" properties, that is, readily 
enter into substitution reactions and very reluctantly 
into addition reactions. As to the compound obtained, it 
behaved quite differently, not entering into substitution 
reactions and easily adding bromine. 

Only 43 years later the veracity of, the German chemist 
was proved. Cyclooctatetraene has no aromatic properties 



* However, the one-electron approximation is inadequate for 
describing the ground state of cyclobutadiene which according to 
more thorough calculations is likely to be singlet. 

§6 



because \ts molecule is not planar but is shaped like 
a tub: 




Moreover, this molecule is believed to have the ability 
of turning "inside out" via a transitional plane structure. 
At the same time it was demonstrated that in a plane 
cyclooctatetraene molecule the bond lengths do not get 
equalized, the double bonds retaining their distinctive 
features. Even if this molecule were plane, it would 
still be nonaromatic. 

In Chapter 4 we shall return to the problem of finding 
out when a system manifests aromatic properties. Just 
now we shall dwell on a very important criterion of 
aromaticity, the Huckel rule. 

The Huckel rule. First, look at Fig. 5 showing the 
graphic method of determining the relative arrangement 
of MOs which was proposed by A. A. Frost and B. Musulin 
in 1953. A regular N-gon is inscribed into a circle so that 
one of its vertices is located at the lowest point. Each 
of the angles of that iV-gon has the corresponding n-MO 
shown on the right side of the figure, the nonbonding MOs 
(e = a) lying on the horizontal line passing through the 
centre of the circle. 

It can be seen in Fig. 5 that monocyclic hydrocarbons 
C n E. n have filled shells (all BMOs are occupied) if the 
number of their it electrons is equal to An + 2, where 
n = 0, 1, 2, .... This is the Huckel rule. It prescribes 
that the aromatic monocyclic compounds should possess 
2, 6, 10 ... n electrons. 

Clearly, neither cyclobutadiene nor (even plane) cyclo- 
octatetraene can be aromatic systems since they possess 
An p n -electrons. 

Non-benzoid aromaticism. From the Huekel theory it 
follows that such systems as cyclopropenyl cation, cyclo- 
pentadienyl anion, cycloheptatrienyl cation (also called 
tropylium cation) and others must exhibit aromatic 
properties. 

The discovery of the two latter ions has an interesting 
story. The C 5 Hj; ion was discovered as early as 1901, but 
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Fig. 5 



only 50 years later didj it attract the attention of re- 
searchers when the first "sandwich" compound, ferrocene*, 



* See p. 28 in the book Symmetry in the World of Molecule^ 
mentioned above, 
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was obtained which proved to have aromatic propertied 
similar to those of benzene. As to the discovery of tropyli- 
um bromide by Merling in 1891, nobody noticed it at 
all, and Merling himself thought it to be a solid isomer 
of dibromocycloheptatriene. At the beginning of 1930s 
Hiickel (not knowing, of course, about Merling's work) 
predicted the existence of tropylium bromide. Only in 
1954 did Dering and Nocks succeed in the intentional 
synthesis of C 7 H 7 Br. 

For the time being we shall drop the subject of aromat- 
ic systems and continue with the analysis of the Hiickel 
method. 

The order of a chemical bond. Many concepts of the 
classical theory of chemical structure were further genera- 
lized and developed in quantum theory of molecules. Here 
we shall touch upon one of such concepts, the bond order. 

Let us leaf back a few pages and look at Eq. (1) which 
provides the mathematical expression of the MO method 
in the LCAO form. When the coefficients c }i are known, 
the order of it bonds can be calculated. The order of the n- 
bond between the atoms / and k is given by the following 
formula: 

N 

Plh= L giCn c hu (13) 

i 

where i is the number of a MO, g t = 0, 1, or 2, the num- 
ber of electrons in the given MO. 

As an example, let us consider the molecules of ethylene 
and butadiene. 

In the former case (see p. 21) 

Pi 2 = 2-l/|/2-l/|/2 = l. 
In the latter case 
p 12 = 2 0.3717 -0.6015 + 2-0.6015 -0.3717 = 0.894 = p M , 
p 23 = 2-0.6015-0.6015 + 2 -0.3717 ( - 0.3717) = 0.447. 

The order of a o*-bond is assumed to be equal to unity. 
Therefore, for ethylene the total order is equal to 

P l2 = 2 
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and for butadiene 

P 12 = -f 34 == 



1.894 and />,,=- 1.447. 



jt-electron density at an atom. In the same manner, know- 
ing the AO coefficients in Eq. (1), one can calculate the 
ji-electron density qj and the charge Qj at the atom /. 
Usually it is done by means of the formulas 



where TV is the number of n-electrons delivered to the 
system by the atom /. 



This way we obtain for butadiene: 
Q l = l _ Ql = 1 _ (2-0.3717 2 + 2- 0.601 5 2 ) = = <> 2 . 



The definitions (13), (14), and (15) of bond orders and 
charges at atoms were proposed by the British chemist 
and physicist C. Coulson. Different definitions of the same 
quantities were suggested at various times by other scien- 
tists (R. Mulliken, K. Ruedenberg, L. Pauling and others). 
We shall not cite or analyse any of those definitions here. 
It is sufficient to note that after some additional sim- 
plifying assumptions all of them reduce to Coulson's 
formulas. 

From the quantities Pjh and q } one can compose a mat- 
rix whose diagonal elements are equal to q } and nondiago 
nal to pj h : 



After Coulson this matrix is called, the matrix of charges 
and bond orders. 

Systems categorized. The electronic structure of hydro- 
carbons is very convenient to describe if the latter are 
subdivided into two classes of alternant and nonalternant 
systems, AS and NAS in the abbreviated form. 



N 



Qj = S Stch, 



(14) 



i=l 



Qj = N-q Jt 



(15) 
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The systems are called alternant if their carbon atoms 
can be conditionally divided into two groups so that each 
atom of one group finds itself directly bonded only to the 
atoms of the other group. Essentially, atoms are subdi- 
vided as follows: in graphic formulas carbon atoms are 
alternately marked with asterisks. If then the marked 
atoms become bonded only to the unmarked ones, the 
system is alternant. For example: 




If such a marking cannot be accomplished, the system 
is nonalternant. For example: 




fulvene azulem 



Note that alternant hydrocarbons have no odd cycles. 

The classification of conjugated systems into alternant 
and nonalternant ones covers compounds both with even 
and with odd number of carbon atoms. However, cer- 
tain properties of even and odd ASs may differ. 

A number of remarkable theorems can be proved for 
ASs; the most significant of them is ... 

... the pairing theorem. Strictly speaking, that is not 
one theorem, but a set of them; below we shall quote a few 
of them. 

1. Orbital energies in even ASs form pairs: + z t and 

— B t . 

2. In odd ASs in addition to paired levels with ener- 
gies ±^i there is also a level with energy equal to a. 

Sometimes the value a is assumed to be a zero reading, 
and the energy is expressed in units of p. Then it can be 
stated that in odd ASs there is a zero level (e = 0). For 
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example, in the benzyl radical (in units of 0):.* 

CH a — 8 1 = e 7 = -2.1010, 
e 2 = e 6 - —1.2593, 
63 = 65= —1, 
e 4 ==e = 0. 

The orbital with e = e is usually called the non- 
bonding MO (NBMO). 

3. If two MOs form a pair and one of these orbitals 
(tyt) is defined by the expression 

itf =2 c Ji<P*-t- -SeSi?*! 
j h 

where the asterisk indicates the AOs (and their coeffici- 
ents) of marked atoms and the small circle the AOs of 
unmarked atoms, the second MO (i|) t r ) can be derived 
from the first one by changing the sign of one of the sums: 

j h 

or \ 

V|5l = — Tj c*i<P? • t- 21 c'hiVk- 

i k 

4. If an AS has a single NBMO, it can be shown that 
the number of atoms in one group (say, the number of 
marked atoms) exceeds by one the number of atoms in 
the other group. Having marked the atoms of the larger 
group by an asterisk, we obtain, for the case of NBMO 

5. The n-electron density q t of a carbon atom in 
neutral ASs, both even and odd, is equal to unity. 

6. In neutral ASs the order of bonds between two 
marked or between two unmarked atoms is equal to zero. 

We have presented only a small fraction of what could 
be told about the structure of unsaturated organic com- 
pounds, even confining ourselves to the Hiickel method. 
Still what has been told clearly shows how useful is this 
simple method of qualitative description of numerous 
organic molecules. During the last three decades the 
HMO started to be interpreted in terms of the graph the- 
ory. This aspect will be covered in the next chapter. 
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In the late fifties the Hiickel method acquired a new 
interpretation on the basis of the graph theory. Since 
that time it became quite popular in chemistry, in partic- 
ular, organic chemistry. Now we shall acquaint the 
reader with the fundamental concepts of the graph theory 
and its relation to the HMO method. 

What Is a Graph? 

Spontaneous definition. In the everyday life we come 
across a graph fairly often, for example, when scrutini- 
zing the map of city streets or highways or railroads, when 
checking up on an electric circuit, when drawing a tech- 
nological flow-sheet, and so on. Here is the most elemen- 
tary definition: a graph is a set of points connected by lines. 

The points are called the vertices, and the lines, the 
edges of a graph. Examples of graphs are given in Fig. 6. 

The degree (or the valence) of a vertex is the number 
of edges converging at that vertex. Thus, the degree of 




(°) (b) ( C ) 

Fig. 6 

vertex A (Fig. 6a) is equal to three and the degree of the 
vertex C in the same figure to two, etc. The vertex of 
degree is called isolated (Fig. 6c, the vertex F), and 
the vertex of degree 1 is called planted or terminal 
(Fig. 6c, the vertex E). 
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When we consider a graph corresponding te'a certain 
real object or process, or just a graph as a mathematical 
object, we take no interest in its metric properties, such 
as the shape of its edge, its length, the exact geometrical 
arrangement of its vertices, etc. What is important to 
know is how many vertices and how many edges a graph 
possesses and how its vertices are connected. If, examining 
two graphs G x and G 2 , we observe that between the sets of 
their vertices there is a one-to-one correspondence, such 
that the number of edges connecting any two vertices of 

C E 
Fig. 7 

G x equals the number of edges connecting the respective 
vertices of G 2 , then we call the graphs isomorphic (see 
Fig. 7). 

Two vertices A and B of the graph G are called adjacent 
if they are interconnected bythe edge AB; in this case 
the vertices A and B are called incident at that edge and 
the edge AB is called incident at the vertices. 

It may so happen that two vertices are interconnected 
by several edges. Then such edges are called multiple 
(Fig. 6c). An edge can connect a vertex with itself; 
such an edge is called a loop (Fig. 66). 

A graph featuring loops and multiple edges is often 
called a pseudograph, and a graph featuring multiple edges 
(but no loops) a multigraph. A graph that possesses 
neither loops nor multiple edges is called simple. Later 
we shall give you examples of some other kinds of graphs. 

More rigorous definition. For those readers who prefer 
mathematical strictness we present a more accurate defini- 
tion of a graph: the finite graph G is the pair [V, E\ 
comprising the following objects: (I) a finite nonvoid set 
V possessing p vertices and (2) a set E possessing q edges, 
i.e. q nonordered pairs of different vertices of V. A graph 
with p vertices and q edges is called a (p, g)-graph. Usual- 
ly a graph is represented descriptively by a diagram which 
is often also called a graph. 
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Subgraph. Certain vertices and edges can be deleted from 
the graph G (Fig. 8). What is left, provided, of course, 
there is something left, is also a graph called a subgraph 



subgeaphs of G' 



Fig. 8 



of the graph G. Here is a more accurate definition: a sub- 
graph of graph G is a graph whose all vertices and edges 
belong to G (Fig. 8). 

Walking around a graph. Let us choose a graph shown in 
Fig. 9 as an example and take a short walk around it trav- 

A * jujj v t 

V 7 Bl € ' 

Fig. 9 



elling along the edges from one vertex to another. The 
initial point of our trip will be the vertex v l and the 
final one v s . 

The alternating sequence of vertices and edges y 1 e 1 y 2 e 2 ... 
...v 7 e 7 v 8 is called a walk, or a (w lt y 8 )-walk. If we take 
some more steps and get back to the initial vertex, our 
walk becomes closed. 

A walk is called a trail if all its edges are different, and 
a path if all its vertices (and edges, naturally) are diffe- 
rent*. A closed trail is called a cycle. The length of a walk 
is the number of edges it contains. In our example the length 
of the (w lt u 8 )-walk is equal to 7. The distance d (v t , v h ) 
between the vertices v t and v h of the graph G is the length 
of the shortest path connecting these vertices. 



* Only the initial and the final vertices can coincide in the 
case of a closed walk. 
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Connectivity of a graph. One of the most important 
properties of a graph is its connectivity. A graph is called 
connected if for any two of its vertices v t and v h there is 
a path from v t to v h , i.e. if we can proceed from v t to 
v k without walking twice along the same edge. 

Forests and trees. The graph theory terminology is not 
yet fully established although the theory is about 250 years 
old. Frequently the terms were taken from the every- 
day life. For example, such words as forest, tree, knot 
have picked up another meaning in this theory. 

To illustrate, a forest is a graph without cycles; if such 
a graph is also a connected one, it is called a tree. Fi- 




Fig. 10 ' 



gure 10 shows four trees which all together form a non- 
connected graph without cycles, i.e. a forest. 

Molecules and graphs. For over 100 years chemists 
depict molecules] by means of so-called structural for- 
mulas, i.e. formulas that do not specify (unless it is done 
deliberately) either bond lengths, or their multiplicity, 
or valence angles, i.e. do not indicate the metric proper- 
ties of a molecule. In terms of mathematics such a for- 
mula is a graph, or a topological graph as they sometimes 
say, in the sense of the word used at the very beginning of 
the book when the topology of a molecule was discussed. 

The only information that structural formulas provide 
is how a given atom is bonded (if at all) to others, or, in 
the words of Butlerov, how "the chemical interaction is 
distributed". Moreover, when chemists began using the 
graph theory, they simplified the representation of a mol- 
ecule as compared to the classical one. In distinction to 
a structural formula a topological graph, as a rule, has no 
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H,0 H 2 C=CH-CH=CH 2 C 6 H 6 

Occasionally the carbon-hydrogen bonds in organic mole- 
cules are neglected: 



As we shall show later, many properties of molecules 
are defined by the shape of their topological graphs regard- 
less of the geometry of molecules or the nature of chem- 
ical bonding in them. 

Graphs, Matrices, Molecules 

Thus, we have established that molecules can be repre- 
sented by means of graphs. To proceed beyond this rath- 
er trivial statement, let us resort again to the mathemat- 
ical graph theory and ask the following question: is it 
possible to represent a graph without resort to a drawing, 
in other words, is it possible to show the relationship of 
the graph's vertices not by means of points and dashes 
but differently, by means of other mathematical objects? 

The answer is yes, and the suitable objects for the 
purpose are matrices. 

Topological matrix of a molecule. There are several ways 
of establishing the correspondence between graphs and 
matrices. We shall dwell only on one of them, the most 
straightforward and convenient. When representing a mol- 
ecule via a graph, we were concerned only with one aspect 
of its structure, namely, we wished to know whether any 
two atoms are chemically bonded or not. Depicting any 
covalent bond by an edge of a graph, we digressed from 




CO 
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such properties of the bond as length, strength, polarity, 
spatial directivity, and multiplicity. 

Let us try to write a topological graph of a molecule 
using a matrix. To do this, we shall first of all enumerate 
the vertices of the graph (i.e. the atoms). Then we shall 
compose a certain matrix A with elements A i} , where i 
and / are the numbers of the vertices (atoms). 

If the atoms i and / are bonded, i.e. the vertices i and / 
are connected by an edge, then A t j — 1; otherwise, A tj = 
— 0. The dimension of the matrix A is equal to the num- 
ber of atoms in the molecule, N . The diagonal elements 



/o l o o o o o o\ 
/ioiooooo\ 



c=c 
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Fig. 11 



A n = 0, i.e. a vertex is not connected with itself and 
the graph has no loops. Such a matrix is customarily 
called an adjacency or a topological matrix: As an example, 
a topological matrix of a styrene molecule is illustrated 
in Fig. 11. 

From vertex to vertex. In the graph theory the matrix 
elements A i} can be interpreted as follows: an A tj is the 
number of unitary walks between the vertices i and /'. 
Then the product of two elements of the matrix A, A Tj A js , 
is equal to 1 if the vertex r is connected with the vertex /, 
and the latter in its turn is connected with s, i.e. between 
r and s there is a walk of length 2 passing through /. If 
there is no such a walk, A T jA js = 0. 

Let us square the matrix A: 

AA = A 2 . 
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Then the element (A 2 ) rs of the matrix A 2 is equal to 

N 

(■^ 2 )rs = £J> A r jAj S . 
j=i 

It is easy to see that the right-hand side of this equality 
represents the number of all walks of length 2 connecting 
the vertex r with the vertex s. Consequently, (A n ) rs is 
the number of all walks of length n lying between r and s. 

First disappointments. Thus we have established a rela-, 
tionship between the graph of a molecule and its topologi- 
cal matrix. Unfortunately, this relationship is not unique 
since the appearance of the matrix A depends on the enu- 
meration mode of the atoms in the molecule. For example, 
the standard enumeration of a graph with three vertices 

1 2 3 

/0 1 

results in the topological matrix I 1 1 

\0 1 

while the enumeration 

1-3 2 

/0 1\ 

originates another matrix I 1 I. 

\1 1 0/ 

In search of invariants. Are there possibly any other 
characteristics of a graph (or its topological matrix) that 
are independent of the vertex enumeration mode? Yes, 
such invariant characteristics do exist. However, they 
can be obtained only after certain refinements of the 
theory. 

Let us consider the following equality: 

Ac ; = XjCj (16) 
(i = 1, 2, . . ., N), 

where c is a column matrix composed of N quantities 
c iii Cj«> • • •> c nii whose quantum-mechanical meaning 
will be clarified later; (x t , 
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of other quantities whose physical meaning is not con- 
sidered so far either. 

The column matrix c ; is called an i-th eigenvector of 
the topological matrix A, and x t an i-th eigenvalue of 
this matrix. The totality of eigenvalues of the matrix A 
is called a spectrum of the graph corresponding to this 
matrix. 

Now let us compose the following matrix from the eigen- 
vectors c t : 

C = (CiCj . . . e N ). 

Surely, the matrix C looks as a row only in the abbre- 
viated notation employed above. Actually, C is a square 
matrix since each of its elements C,- is a column matrix. 
Using the matrix C, Eq. (16) can be written as follows: 

AC = CX, 
where X is the diagonal matrix 




Bounds of the spectrum of a graph. The spectrum of 
a graph possesses a remarkable property: it is independent 
of the way the vertices are enumerated. Besides, the spectrum 
has still another important property: if D max is the 
maximum degree of a graph's vertex, then 

This result is essentially the same as the theorem de- 
rived by F. Frobenius in the 19th century. 

The quantities x t may or may not reach their extreme 
values, but they never fall outside the interval [— D max , 

As a rule, x t are not integers (in contrast to D max ), but 
some graphs have spectra consisting of whole numbers 
(see Fig. 12, the values of x t are given in braces). 

In conjugated molecules D max = 3. Therefore all 
Hiickel energy levels lie in the interval [—3, +3]. Simi- 
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larly, in 4he linear polyenes and annulenes D mnx — 2 and 
so the extreme values of these levels are — 2 and +2 in 
these molecules. 



{t.-i} 

Dmax = l 



■A, 
{2,-1,-.} 



{2,0,0,-2} 




{2,1, 1,-1.-1, -2} 




{2,1,0, 0,-1,-2} 




{2, 1, 1, 0.-1,-1 - 2 } 



Fig. 12 



Isospectral graphs. When discussing the spectrum of 
a graph we stressed its invariance relative to the vertex 
enumeration mode. Nevertheless some topologically non- 
equivalent graphs can sometimes possess an identical 
spectrum; such graphs are called isospectral, or, in math- 
ematical literature, cospectral. The spectra of topological 
graphs of 1,4-divinyl benzene and 2-phenyl butadiene 
molecules (Fig. 13) can serve as examples. 

How to determine the spectrum of a graph? Let us rewrite 
Eq. (16) in the form 

(A — Xjl)c,- = 0, (17) 

where I is a unitary matrix, and is a zero matrix. 

The attentive reader has guessed right: we are nearing 
a secular equation. Indeed, Eq. (17) has nontrivial solu- 
tions if the following condition is satisfied: 



det | A — x t l | - 0. 



This secular determinant is nothing but a polynomial of 
x t which is characteristic of a given graph G. This poly- 
nomial is appropriately called a characteristic polynomial 

molecules 




their graphs 





graph spectrum 

{ 2.214; 1.675; 1.0; 1.0; 0.538;"* 
-0.538; -1.0; -1.0; -1.675; -2.214 } 



Fig. 13 



of the graph G and is denoted by the symbol P G (x): 
P G (x) = det \A — Xil\=a x N +aiX N ~ l + . ..-\-a N = 

N 

= 2 a n x"-". (19) 

Tl=0 



Therefore, Eq. (18) can be replaced by the following one: 
P G (x)= 2 a n x»-» = 0. 

n=0 

Thus we have derived an algebraic equation of N-th 
degree, which has, in accordance with the fundamental 
theorem of algebra, N roots. However, sometimes it is 
not as important to calculate the roots as to analyse 
their properties, e.g. to find out how many positive, 
negative, zero roots the equation has, and so on. In such 
cases we have to find the coefficients of the characteristic 
polynomial. But how? 
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Sachs' theorem. In 1962-64 the German mathematician 
Horst Sachs (GDR) proposed an imaginative method of 
determining the coefficients of characteristic poly- 
nomial*. Here is the description of his method. Suppose a 
molecule can be represented by the following graph G: 




Let us "expand" this graph into subgraphs whose ele- 
ments are either cycles or fragments of the type o — ° 

Such subgraphs are called the Sachs graphs (G s ). 

In our case we can distinguish in the given graph, for 
example, 5 Sachs graphs shown in Fig. 14. 

To proceed, it is essential to know: 

(a) how many components in all the given Sachs graph 
c (s) has and 

(b) how many cyclic components r (s) it features. 
The numbers c (s) and r (s) are easy to determine. For 

example, for the graph G s „ c (s) 
Fig. 14). 

A / A 
V V 

Fig. 14 

Now we can formulate the Sachs theorem: 

The coefficients of a characteristic polynomial of a certain 
graph G can be obtained from the following formulas: 

a =i, 

a„= J\ (-l) c < s >-2 r<!, > (1 (20) 

* The first to suggest a topological method lor constructing 
a characteristic polynomial was C. Coulson (1950). 



= 2 and r (s) = 1 (see 
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where S n is the set of all Sachs graphs with n vertices and 
the summation is carried out over all Sachs graphs with 
n vertices. If S n is a void set, then a n = 0. 

Let us illustrate this theorem by specific examples. We 
shall begin with a simple "diatomic" graph (N = 2) of 
the type 



C 

This graph is a Sachs graph and it cannot be splitted 
into other Sachs graphs. 

Clearly, S t is a void set since there are no Sachs graphs 
with one vertex, and therefore a x = 0. Note that iSx is 
always void in accordance with the definition of the 
Sachs graph and % = all the time. v 

Next, iS 2 = jo — oj , i.e. there is one Sachs graph 

with two vertices; it consists of one component and 
has no cycles. Therefore, c (s) = 1, r (s) =0, and a 2 = 
= (— 1) 1 2°= — 1. Consequently, the characteristic poly- 
nomial takes the form (see Eqs. (19) and (20)). 

P G (x) = aoX 2 + a^x + a % = x 2 — 1 . 
Another example (N = 4): 




This graph corresponds to the following Sachs graphs 
with indicated indices c (s) and r (s): 



MI. / c(s) "'- r(s) "° 
MA} .w-<. -(•)-« 

c(s) =2 , r(s) =0 
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Then in accordance with the Sachs theorem: 
a 2 = (-l) 1 2° + (-l) 1 2 +(-l) 1 2°-f(-l) 1 2 = -4, 
a 3 = (-l)^ 1 = -2, 
a 4 = (_i)i 2 o = 1, 
p G ( x ) = x 4 — 4x 2 — 2* + l. 

Of course, the Sachs method becomes practically in- 
convenient for large molecules whose graphs are unwieldy, 
and therefore scientists proposed other methods of deter- 
mining the coefficients of characteristic polynomial. 
Still, the Sachs method proved to be quite popular. 

Corollary of the Sachs theorem. A few interesting con- 
clusions can be deduced from the Sachs theorem. We have 
already obtained corollary one:* 

a x = 0. 

Corollary two: the sum of the roots of a characteristic 
polynomial, that is, the sum of the eigenvalues of the 
topological matrix of a molecule, is equal to zero:?* 

Xi + + • • • +x N = 0. 

Corollary three: the sum of the squares of roots of 
a characteristic polynomial is equal to ( — 2) a 2 : 

x\ + x\ -j- . . . + x% = — 2a 2 . 

Corollary four: the sum of the cubes of roots of a charac- 
teristic polynomial is equal to (— 3)a 3 : 

x\ + x 3 t + . . . +x%= -3a 3 . 

Corollary five: the number of Sachs graphs of the type 

« — o is equal to the number v of the edges in the 

initial graph G, whence it follows that a 2 = — v. 

Corollary six: a 3 = — 2n 3 , where n 3 is the number of 
three-membered cycles in the initial graph G. 



* In the case of compounds with heteroatoms 5 X becomes a 
nonvoid set and a, ^= 0. 

** The sum of the eigenvalues of a matrix whose all diagonal 
elements are zeros is always equal to zero. 
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All these corollaries are not difficult to prove, and we 
l#%ve the proofs for our readers. Besides, new corollaries 
"Can be derived from them, for example, these two equali- 
ties: 

x\ + x\+ ...+x% = 2v, 
x{ + x\+ ... + x a N = Gn 3 . 

The number of corollaries can be increased for during 
the centuries of studies of algebraic equations mathemati- 
cians proved a great number of theorems aiding to find 
out the properties of their roots without solving the 
equations themselves, knowing only their coefficients. 

But let us leave the graph theory and algebra for a while 
and return to chemistry. We do this to show the connection 
between the material discussed above and the molecular 
structure theory. 



Two Approaches to One Problem 

Topological orbitals... Probably the smartest of our 
readers have already appreciated the close relationship 
between the graph theory and the Hiickel method. Actual- 
ly in quantum chemistry there is a great variety of prob- 
lems in which the Hamiltonian of a molecule can be 
written in a matrix form as a one-valued function of the 
topological matrix of that molecule: 

H = H (A). 

In the simplest case (the Hiickel method) this relation 
takes the form 

H = oI + pA, (21) 

where a is the Coulomb integral and p is the resonance 
integral. 

With Eq. (21) in mind, compare the Hiickel equa- 
tion (9) with Eq. (16) defining the spectrum of topological 
graph of molecule: 

HCj - SCjE; 

and 

Acj = CjXj. 
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First »f all it should be noted that for some types of 
the matrix S (see below) the eigenvectors c ; of the topo- 
logical matrix of the molecule coincide precisely with those 
of the Hiickel method, i.e. with the series of coefficients 
defining the molecular orbitals calculated by this method. 

The coincidence of the eigenvectors of the topological 
matrix and Huckel's MOs follows from the fact that H 
and A commute, i.e. [H, A] = 0. 

That is why the molecular orbitals calculated by the 
Hiickel method are called topological. 

... and their energy. But this is not enough. Let us write 
the overlap matrix S in the form 

S = l4-oA, (22) 

where a is the overlap integral. 

This notation corresponds to the following assumption: 
if the atoms are chemically connected by a "direct" cova- 
lent bonding, their orbitals overlap and the values of the 
overlap integrals are always equal to the same magni- 
tude a. If the atoms are not bonded directly, then their 
orbitals are supposed to be not overlapping and the 
respective elements of the matrix S to be equal to zero. 
It follows from the norming condition that the diagonal 
matrix elements are equal to unity. Therefore, the ma- 
trix S expressed in the form of (22) is obviously topological. 

Now, substituting expression (21) for H and (22) for S 
into the Hiickel equation (9), we obtain 

(H — eiS)c,- = (al + 0A — e,I — e^A)^ = 0, 

or, taking into account the properties of a unitary matrix 
Acj = CjX; and Ic; = Cjl = Cj, 

(a -f- $x t — e,- — 8jax,)cj = 0. 

Whence, equating the expression in parentheses to 
zero, we get: 

Thus, we have expressed the orbital energies of the 
simple MO method in terms of the eigenvalues of the topo- 
logical matrix. If we now pass over to the simplest ver- 
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sion of the Hiickel method, that is, apply the«ZDO ap- 
proximation, we have 

e; = a]+ P-T;. 
If we also recall that in this simplest version of the 
HMO method the value of a is often taken to be the zero 
reading and p a unit of energy, then 

We have thereby clarified the meaning of eigenvalues of 
a topological matrix of a molecule: these values, when 
expressed in units of p, coincide with the- orbital energies 
of the HMO method. In other words, the Hiickel [method 
determines the MO energy values specified only by the 
molecular topology. 

It must be emphasized, however, that this coincidence 
is due to the peculiar nature of the Hiickel Hamiltonian, 
with the short-range forces being dominant in the effective 
potential. 

Usually the Hiickel topological MOs are classified into 
bonding, antibonding, and nonbonding ones, depending 
on the sign of the respective orbital energies. The MOs 
corresponding to the values x t > are called bonding, 
the MOs for which x t < antibonding and those with 
x t = nonbonding. 

The numbers of bonding, nonbonding, and antibonding 
orbitals are denoted by N+, N , and N. respectively, 
and are related to the total number N of atoms in a given 
graph (a conjugated molecule): 

N + + N + N. = N. 

The quantity (N+-N .) is called the graph (molecular) 
signature. 

The graph-theoretical methods are extensively used in 
the theory of conjugated systems. It is only natural since 
in this area of chemistry the HMO method was used most 
frequently. In the next chapters we shall dwell on some 
interesting results derived in the topological theory of 
conjugated and aromatic systems. 



Chapter Three Graphs 

of Conjugated 
Molecules 



In this chapter we shall discuss some significant results 
obtained in the theory of conjugated organic molecules by 
the graph-theoretical methods. 

Coloured Graphs 

How to colour a graph. In Chapter 1 we subdivided the 
conjugated systems into two classes, alternant and non- 
alternant. It is obvious that such a subdivision is based 
on topological features of molecules. Indeed, classifying 
a system as alternant, we digress from its metric proper- 
ties and emphasize the way in which the atoms are inter- 
related. In one case the atoms are connected so that they 
can be categorized in the manner described on page 31, 
in another case their connection is quite different. How to 
formulate such a classification in terms of the graph theo- 
ry? To do this we have to colour the graph vertices so that 
identically coloured vertices are not connected by any 
edges. It should be pointed out that in the general case 
the problem of colouring a graph with an arbitrary num- 
ber of colours, so that the adjacent vertices are coloured 
differently, is extremely difficult. The solution is found, 
however, for some specific cases. We shall be concerned 
with the problem of two colours*. 

The graph that can be marked in two colours is called 
bipartite, or bichromatic. Hereafter we shall designate 
the vertices of a topological molecular graph marked in 
one colour by an asterisk, and in another colour by a small 
circle, i.e. we shall keep to the designations adopted in 
Chapter 1. Figure 15 gives examples of bipartite (a) and 

* For all that, in principle, all chemical graphs of organic 
molecules can be coloured in three colours. 
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nonbipartrte (b) graphs. The former correspond to alter- 
nant, and the latter to nonalternant systems. 

Biparticity criterion. Look carefully at Fig. 15 to see 
that the bipartite graphs shown in it break up into the 
subgraphs possessing only cycles with even number of 




Fig. 15 



vertices or not possessing cycles at all. This observation 
can be generalized and the following theorem proved: 

The graph is bipartite if, and only if, it has no odd- 
membered cycles (rings) among Us subgraphs. 

The formal proof of this theorem was given by D. Konig 
in 1936. 

Topological matrix of a bipartite graph. Bipartite graphs, 
just as the corresponding alternant systems, possess a num- 
ber of remarkable properties. In particular, their verti- 
ces can always be enumerated so that the topological 
matrix is simplified and reduced to the block form 

/0 B\ 

a = (bo)- < 23 > 

where B is a certain matrix which is a block of the matrix 
A, and Bis a transpose of the matrix B, i.e. the matrix 
obtained from the matrix B by interchanging its rows and 
columns: B ik -*■ B hi . For example, the topological mat- 
rix of a benzene molecule whose atoms are enumerated 
as shown (first the marked atoms are enumerated and 
then the remaining ones) 



1 




5 
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takes the .form! 



A = 



Oil 1 





1 1 





1 1 


1 1 





1 1 





1 1 






Pairing theorem again. And how is the characteristic 
polynomial of a bipartite graph look like? Recalling the 
Sachs theorem one can easily write the expression for 

Pg (*)= 

P G (x) = x" + a 2 x N - 2 + a k x N - i + .. . . 

Indeed, since a bipartite graph has no subgraphs with 
odd-membered cycles, the set of the Sachs graphs S n , 
where n is odd, is void (S n = 0) and a n = 0. 

Therefore, when N is even, 

P G (-x) = P G (x), 
and when N is odd, 

P G (-*) = ~P G (*)• 

But in either case, if x t is a root of the characteristic 
polynomial 

P G (**) = 0, 

then (— x t ) is also its root. 

On the top of that, there is at least one zero root 
(x t = 0) for odd N (odd ASs). 

Now compare the conclusions deduced to what was said 
in Chapter 2 and you will see that what has just been pre- 
sented is nothing else but the pairing theorem formulated 
in terms of the graph theory. 

Proceeding with our arguments, we can demonstrate 
that the eigenvectors of the topological matrix of a mole- 



cule which correspond to the bonding MOs have the fond 



P BMO_ 



c*u 

ch 



c*. 



c P+qi . 

(the vertices 1, 2, . . ., p are marked, and the vertices 
p + 1, . . ., p + q = N are not marked). 
For antibonding MOs: 



tABMO _ 



r* 

c Ph 



c P+gi 

and for nonbonding MOs (see item 4 of the pairing theo- 
rem described earlier): 



pNBMO _ 



Cpl 







-0- 



(24) 



Thus, the complete matrix composed of AO coefficients 
(if NBMOs are absent) tak.es the following form in the 
case of ASs: 



C = ( C BMO. C ABM0) ; 



(v -v)' 



where U and V are the blocks of the matrix C. 
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The matrix composed of the eigenvalues of topological 
matrix of a bipartite graph whose vertices are properly 
enumerated also takes the sufficiently simple block form: 




where the matrix Y = 2UBV = 2VBU is a diagonal 
matrix. 

Bond orders and topology of a molecule. In 1955 the 
British physicist G. Hall and in 1958 the American scien- 
tists N.S. Ham and K. Ruedenberg found that the bond 
orders, whatever was the method of their definition 
(according to Goulson, Mulliken, Pauling, or otherwise), 
can be expressed via the topological matrix of a molecule. 
Thus, the matrix of bond orders and charges, constructed 
according to Coulson is 



As has already been said, the system representing an 
even alternant hydrocarbon has a topological matrix of 
the type (23). Then, in accordance with the Hall theorem 
the matrix of orders of bonds between marked and un- 
marked atoms (P*°) is expressed by the following topolo- 
gical formula: 



Without going into particulars over these formulas we 
shall formulate the general statement: 

The bond order depends on the topology of a mole- 
cule, and since the bond order correlates with the bond 
length, the geometry of the molecule is to some 
extent defined by its topology. 

Let Us Count Orbitals! 

The foremost problem of theoretical chemistry is the 
stability of a chemical compound. TheJ problem is ex- 
tremely difficult because it involves two aspects, thermody- 
namical and kinetic. Nevertheless, in some special cases 
a number of rules can be established from which predic- 




P*° = (BB) 1 / 2 B. 
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tions can be made. For example, stability (or instability) 
of some n -electron systems is often determined to a large 
extent by their topology. 

We shall begin with the problem that is essential for 
chemistry, the determination of the number of bonding, 
antibonding, and nonbonding MOs in a certain jx-electron 
system. From the standpoint of the graph theory the prob- 
lem can be stated differently: how many positive, nega- 
tive, and zero eigenvalues are there in the spectrum of 
the graph mapping the topology of the given n-electron 
system? If the number of n centres in such a molecule is 
equal to N (a graph with N vertices), the total number of 
n-MOs of all types is N. Consequently, if there are no 
nonbonding MOs (we shall explain later on how to estab- 
lish their] presence), the task becomes easier: we have to 
determine only the number of, say, antibonding MOs. 
The latter can be done by means of the graph theory, 
using a comparatively simple method. 

Annihilate a graph! Suppose we consider a certain graph 
G with N vertices «ach of which is denoted by the letter 
Ki, with the index i labelling the vertices. Now let us 
delete one of the vertices together with the edges incident 
at it. This way we get a new graph that has one vertex 
less: 

Repeating this operation, we shall have the graph 
G N . 2 . Suppose we continue deleting one vertex after anoth- 
er until there will be nothing left of|the graph, or, as the 
mathematicians say, we have the void graph G . The 
sequence of graphs obtained in successive deleting of 
vertices from a certain initial graph, can be represented 
as follows: 

{G N , G N _ t , . . . , GiG }, 

where G N -! = G N —K N , G N . 2 = G,,.,-^-!, etc. 

This ordered set of subgraphs of a graph G is called the 
expansion of G, or the decomposition-sequence of G. 

Some more designations. Making notations is an indis- 
pensable — even though boring — job and we shall occupy 
ourselves with it now. First, we shall denote the number 
of positive, zero, and negative eigenvalues in the spectrum 
of the graph G n by N+ (n), iV ( re )> an( l N - ( n ) respective- 
ly. Then, the symbol D N . n will be used to denote the 
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determinant of a topological matrix of the graph G N . n 
(n = 0, 1*. . ., N). And finally, the sequence of numbers 
{D N ; D N . . ., D } will be designated by the sym- 
bol W N ]. 

Number of ABMOs. Now we can formulate the theorem 
on the number of ABMOs: 

The number of ABMOs (N _ (n)) is equal to the number of 
sign changes in the sequence W N ]. 

How many ABMOs are there in pentalene? It is common 
knowledge that the best way to clarify the meaning of 
a theorem (especially this one, quoted without proof and, 
I regret, with some simplifications) is to give an exam- 
ple. As such an example, let us consider a pentalene 
molecule: 




First, let us annihilate the graph of this molecule 
(Fig. 16) and then compose a topological matrix for each 

<o-Ov-<x- 

«. O, G, 

— — ^ - ^ 

G - G - G, G 2 G, G„ 

Fig. 16 

of the subgraphs thus obtained and calculate its determi- 
nant. We advise the reader to do it for himself as a useful 
exercise* . 

Here are the values of the determinants: 

[D N ] = (0, -2, -1, | 2, 1, | -1, 0, | 1). 

The vertical lines indicate the places in the sequence 
at which the sign changes. Tbere are three such places, 
so W.(8) = 3. 



f for a void graph D„ .i X. 



Some qualitative conclusions can easily be cbrawn from 
the results obtained, in particular, the tendency of the 
molecule to attract electrons to the NBMO. Indeed, 
a pentalene dianion is fairly stable. 

ABMOs of benzene. Let us turn now to a benzene mole- 
cule as an example ofj an alternant hydrocarbon. In this 
case the expansion of the initial graph has the form illus- 
trated in Fig. 17. 

0-0-n-v*— - 

G e G s G 4 G 3 G 2 G, G 

Fig. 17 



Reasoning as above, we obtain the following sequence 
of values of the corresponding determinants: 

[D N ] = (-4, 0, 1, 0, -1, 0, 1). 

Whence N -(6) = 3, which agrees with the results ob- 
tained above. 

Thus, the described method is applicable both to ASs 
and to NASs. 

Are there NBMOs in a molecule? In some cases the in- 
formation on the stability and reactivity of a hydrocarbon 
can be assessed by determining whether the compound has 
singly populated NBMOs. In terms of the graph theory the 
presence of NBMOs is determined by establishing the 
necessary and sufficient conditions for the existence of 
zeros in the graph's spectrum. Below we cite a number of 
important theorems allowing the presence of NBMOs in 
a given AS to be established. 

A few words on N+. The quantity N + , i.e. the number 
of bonding orbitals, can be found from the characteristic 
polynomial of a molecule in the following way: 

a N ~ a N-l = • • • = fljv-JVo + 1 = 0» 
a N -No =7^ 0, 

Ch (a N , a N - x , . . ., a u a ) = N+, 

where Ch denotes the number of changes of sign in the 
given sequence of coefficients. As an example, consider 
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fulvene*(JV = 6) (see p. 31). Its characteristic polynomial 
is 

P (x) = x 6 - 6x 4 + 8x 2 - 2x - 1, 

or 

x 6 + 0-x* — 6-x i + 0-xZ + 8x z — 2z — i. 

Then 

Cft (-1, -2, +8, 0, -6, 0, +1) = N+ = 3. 

In addition, 7V = and N. = N — N — N+ = 6 — 
_0 -3 = 3. 

Six Remarkable Theorems 

Theorem one — on the number of zeros. Suppose a certain 
alternant hydrocarbon possesses p marked and q unmarked 
atoms. Then the number of zeros (N ) in the spectrum of 
the topological graph of this molecule obeys the in- 
equality 

No>P — g- 

The following corollary can be derived from this theo- 
rem: if an alternant hydrocarbon possesses N carbon 
atoms (N n-centres) and there are no 4m-membered 
cycles, then 

N = N - 2t, 

where t is the maximum number of double bonds that can 
be written in the structural formula of the given hydro- 
carbon. 

Let us consider benzene, for example. Its molecule has 
no 4m-membered cycles, TV = 6, t = 3, whence ./V = 0. 

Generally speaking, when a conjugated^system possess- 
ing no 4/n-membered cycles can be expressed graphical- 
ly by the Kekule structure, i.e. the structure with alter- 
nating double and single bonds, for which N — 2t all the 
time, the system has no NBMOs which contributes to its 
stability. But if this compound cannot be written in the 
Kekule form, then it is a radical (Dewar and Longuei- 
Higgins rule). 

Theorem two — on deletion of two vertices. If in the 
topological graph of an alternant hydrocarbon there is 



57 



a vertex of order one, we can delete it together r with the 
neighbouring vertex and the edges incident at the latter 
without changing the number of zeros in the graph's spec- 
trum (see Fig. 18a, b). 

Theorem three — on deletion of four vertices. Briefly, 
this theorem can be formulated as follows: if a trail con- 




Fig. 18 



sisting of four vertices is replaced by one edge, the num 
ber of zeros in the graph's spectrum remains the same 
(Fig. 18c). 

Theorem four — on destruction of a peripheral ring. If 

the topological graph of an alternant hydrocarbon has 
a peripheral four-membered ping, two vertices and four 
edges can be deleted from that ring without changing the 
number of zeros in the graph's spectrum (Fig. 18a!, e). 

Theorem five — almost obvious. If the topological graph 
G of a molecule can be represented in the form of two siib- 



graphs Connected as follows: 



and N (d) = 0, then N (G) = N (G a ) (Fig. 18/). 

This seemingly trivial theorem has a practical signific- 
ance: the introduction of the vinyl, phenyl, and other 
groups for which N = can neither essentially stabilize 
an unstable molecule, nor essentially destabilize a stable 
one. 

Theorem six — on one important fragment. If the topolo- 
gical graph of a molecule has the fragment 




then JV >0 (Fig. 18/). 

More about a NBMO. Now let us get back to the vector 
(24) composed of the AO coefficients appearing in the ex- 
pression for NBMO. 

Since in the case of NBMO the values of x t are equal to 
zero, for ASs we obtain from Eqs. (16), (23), and (24) 

Be* = 0. (25) 

In an alternant hydrocarbon each unmarked atom can 
be connected at most to three marked ones: 

/C?(e?o) 

(cfo)Cf-q/ 

The AO coefficients appearing in the expressions for 
the given NBMO (p are given in parentheses. 
Then, from Eq. (25) we get 



cfo -r c% + Cho = 0. 



(26) 



As an example, let us consider the allyl radicaj \GR = 
= CH - CH 2 ) whose topological graph has the form 



Let us enumerate its vertices so that the number of 
labelled vertices would exceed the number of unlabelled 
ones: 



* o * 

1 3 2 



The topological matrix of this alternant system is 
then written in a block form as follows; 

'0 0il N 
1 1 

JTio, 

i.e. B = (J) and B = (1 1). 

The matrix of the coefficients by the marked atoms for 
a NBMO of the allyl radical 



r>HBMO 



gives the following expression when multiplied bv B 
(see Eq. (25)): F 



(1 1) I ■ = C *0 + C 20-=0. 



In the general case Eq. (26) is valid fo r each unmarked 
atom of odd AS. 

Eq. (26) makes it possible to find the AO coefficients 
of marked atoms (the remaining coefficients being equal 
to zero) in the expression for NBMO. Thus for the allyl 
radical examined above we obtain: 



"Redundant" Symmetry 

Graph and symmetry of a molecule. The symmetry of 
a molecule is not manifested explicitly in its graph. For 
example, both graphs shown below correspond to a ben- 
zene molecule: 




Nevertheless, we can introduce the concept of a ... 
... group of a graph. Let us consider a noncyclic graph 
comprising three vertices: 



First, we enumerate the vertices of the graph. It can 
be done in several (more precisely, in 3! = 6) ways, two 
of which are depicted below: 

/ 2 3 I 3 2 

Each vertex enumeration mode corresponds to a definite 
permutation and its own topological matrix of the graph. 
Thus, graph G x is associated with the topological matrix 

/0 1 ON 

A,= 1 1 

\0 1 Oy 

and graph G 2 with 

/0 IN 
A 2 = 1 
\1 1 Oj 

The vertex permutation can be conditionally written 
as follows (for graph G 2 ): 

Sometimes the upper row is omitted. 

Such a notation indicates that vertex 2 in graph G ly 
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whose vertex enumeration is adopted as standard, is 
replaced by vertex 3 in graph G 2 , and vertex 3 in graph G 1 
by vertex 2 in graph G 2 . Vertex 1 is not involved in per- 
mutations. 

The set of permutations of N vertices, or more general- 
ly, N objects, forms a group called a permutation or 
symmetrical group (S N )*. 

In so far as there is a group we can also speak of its 
representation, i.e. of a group of square matrices obeying 
the same group laws as the permutations themselves. 

How to compose a matrix representing a given permuta- 
tion? To answer this question, we shall return to our 
example. The permutation matrix, denoted by the sym- 
bol r(/ > 1 ), is composed as follows: the matrix element lo- 
cated at the intersection of the second row and the third 
column is assumed to be equal to unity since vertex 2 is 
replaced by vertex 3. Accordingly, the element located at 
the intersection of the third row and the second column 
is also equal to unity for vertex 3 replaces vertex 2. 
After the permutation vertex 1 retains its position, and 
therefore the element r u (/ > i) of the matrix r(Pj) is also 
equal to unity. The remaining elements of that matrix 
are assumed to be equal to zero. Thus, the sought matrix 
T(Pi) corresponding to the permutation P 1 has the form: 

/l 0\ 
T(P l ) = 1 
\0 1 Oj 

Similarly, the matrix 

/0 1 N 
r(/> 2 )= 1 

Vl 0/ 

corresponds to the permutation 

1 2 3^ 

P 2 = \ 



/12 8X 
V3 2 i)- 



* The readers are supposed to be acquainted with the founda- 
tions of the group theory. The most general concepts of this 
theory are presented in our booklet Symmetry in the World of 
Molecules, Mir Publishers, 1979. The permutation groups are 
well described in the excellent popular book by I. Grossman 
and W. Magnus Groups and their Graph?, Random House, 1964. 
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The fallowing transformation of a topological matrix 
of a graph plays an important part in the graph theory: 

A' = f (P)Ar (P), (27) 

where f (P) is the transposed matrix. 

For the first considered permutation of the graph's 
vertices the transformation (27) results in a new topo- 
logical matrix differing from A: 

A'-f (P,)Ar (p,) = 

(\ ON /0 1 ON /l 0\ /0 IN 
= 001101001^001 |^A 
\0 1 0/ \0 1 0/ \0 10/ \1 1 0/ 

However, the same transformation performed on the 
matrix V(P 2 ) leaves the topological matrix A invariant: 

/0 IN /0 1 ON /0 IN /0 1 ON 
A' = ( 1 0(1 1 ( 1 = 1 1 U A. 
\1 0/ \0 1 0/ \i 0/ \0 1 0/ 



Thus, in the latter case the permutation of vertices P 2 
leaves the topological matrix A invariant with respect 
to the transformation (27): 

A = f (P) Ar (P). (28) 

The latter transformation is called the automorphism 
of a graph. The totality of all automorphisms of a graph 
forms a group (a group of graph's automorphisms) which 
is in fact a subgroup of its vertex permutation group. As 
we have seen already in the considered example, far from 
all permutations enter a group of graph's automorphisms. 

More rigorous approach. Now let us formulate the ob- 
tained results more strictly. In Chapter 2 we introduced the 
concept of isomorphism of a graph. The isomorphic map- 
ping of one non-oriented graph onto another is a one-to- 
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on mapping of vertices and edges of one graph onto 
the respective vertices and edges of another, so that the 
incidence relations are maintained. Two graphs are iso- 
morphic if there is an isomorphic mapping of one graph 
onto the other (see, e.g., Fig. 7). 

The automorphism of a graph is the isomorphic mapping 
of that graph on itself. For example: 




Gj G 2 



The totality of all automorphisms of a graph forms 
a group. This group is a group of vertex permutations and 
is called a vertex group of a graph, or just a group of a 
graph. When edges are interchanged instead of vertices, we 
deal with an edge group of a graph. When graphs repre- 
sent molecules, both (vertex and edge) groups are isomor- 
phic provided we do not use graphs with loops and mul- 
tiple edges. 

Additional symmetry. When certain organic molecules 
were studied by the Hiickel method an interesting and 
significant property of their topological Hiickel graphs 
was found: the total group of such a graph can consist of 
a greater number of elements than a conventional geometrical 
group of symmetry. 

Let us consider, as an example, a diphenyl molecule 
which can be represented by the following graph: 
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When tl\e vertices of the given graph are enumerated as 
shown, its topological matrix takes the form 
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The geometrical symmetry group of a topological graph 
of a diphenyl molecule includes four elements: the identi- 
ty transformation E and three symmetry axes of second 
order C"*\ C?\ and C' z >. 

We omit the plane of the drawing as a symmetry ele- 
ment because the reflection in this plane produces no ver- 
tex permutations, reducing thus to the element E. The 
four symmetry elements of the graph which were men- 
tioned above form group D 2 and correspond to the follow- 
ing permutations: 



^(C ( e z) ) = ( 



1 2 3 4 5 6 7 8 9 10 11 12> 
2 3 4 5 6 7 8 9 10 11 12, 
1 2 3 4 5 6 7 8 9 10 11 12> 
3 4 1 2 7 8 5 6 9 10 11 12, 
2 3 4 5 6 7 8 9 10 11 12> 
43218765 10 9 12 
2 3 4 5 6 7 8 9 10 11 12> 
436587 10 9 12 



;)• 

12\ 
11 J' 

m 

iij' 



We advise the reader to construct the matrices of these 
permutations for themselves. 

Now let us turn one (for example, lower) ring of a di- 
phenyl molecule through the angle of 180° keeping the 
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other ring fixed, which corresponds to the permutation 
/l 2 3 4 5 (5 7 8 9 10 11 12\ 
\\ 43258769 10 11 12 J " 

The matrix of this permutation satisfies Eq. (28), i. e. 
the total group of the graph turns out to be wider than 
group D 2 . The new symmetry element (k) of the graph to- 
gether with the four others mentioned above form a group 
which is isomorphic to group D id . This fact is essential 
for understanding the features of the energy distribution 
of a diphenyl molecule calculated by the Hiickel method. 
If we adopt the geometrical group D 2 of diphenyl as its 
group of Hiickel graph and assign an irreducible represen- 
tation (IR) of this group to each orbital energy level, we 
shall obtain the following table of orbital energies 
(in P units): 

MO energy IR ( f C 2 group IR of Did group 

—2.278— B l B 2 

— 1.891 — A t A 1 

—1.317— B 1 B t 

—1.000 B 2 , B 3 E 

—0.705— A! A t 

0.000 

+0.705— B l B 2 

-f 1.000 B 2 , B 3 E 

+ 1.317— A 1 A t 

+ 1.891— B l B 2 

+ 2.278— A t Ai 

The levels B 2 and B 3 turn out to be degenerate. The 
origin of this degeneration is not clear if one proceeds only 
from group D 2 . However, if one proceeds from the broad- 
ened group D id of the graph, the degeneration becomes 
clear: this group contains a two-dimensional IR (E), ac- 
cording to which certain MOs of diphenyl are transformed. 

Similar results were also obtained for some other mole- 
cules. For example, the Hiickel graph of a triphenylmeth- 
yl radical has a geometrical symmetry group D 3 (if, 
again, the reflection in the plane oi the molecule is omit- 
ted), whereas its total group is isomorphic to the sym- 
metry group of octahedron (O h ), i.e. belongs to the class 
of cubic groups! Initially this result was obtained by 
means of a computer, but later the problem of additional 
symmetry of a topological graph of a molecule was consid- 
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ered in *njore general terms by V. Ya. Bespalov, research 
worker of Leningrad University. 

The importance of "through-space" and "through- 
bond" interactions. What is the cause of additional sym- 
metry? The fact is that the Huckel method does not take 
into account "through-space" and "through-bond" interac- 
tions. We shall clarify this using a diphenyl molecule 
again as an example. Let us take into account some 
interactions between the atoms that are not directly 
related (shown by dotted lines): 



This graph will be associated with group D 2 , but with 
different topological matrix. 

We would point out once more that the presence of addi- 
tional symmetry is associated with the specificity of the 
Huckel method, namely, with the fact that "through- 
space" 4 and "through-bond" interactions are neglected. 
Incidentally, a similar "mysterious" degeneration (or 
quasidegeneration) can be observed in some other methods 
due to the same reasons. 



Chapter Four Topology 

of Polyene 
Ribbons 



In this chapter we shall resume discussing the struc- 
ture of aromatic compounds. Admittedly, our narration 
cannot lay claim to any comprehensive analysis of the 
problem, which is feasible only in a special monograph. 
Our objective is more modest. We are going to discuss the 
basic aspects and difficulties of the modern theories of 
aromaticity and to show to what extent the concepts of 
chemical topology can be employed in investigating cy- 
clic conjugated systems. 



Fate of a Notion 

In search of aromatic annulenes. The detailed analysis 
shows that the notion of aromaticity is far from being 
simple and unambiguous as it may seem at first glance. 
We shall list here some complications. 

When at the beginning of 1930s E. Hiickel formulated 
his famed (An + 2) rule, only two examples were known 
which confirmed it, namely, benzene and a C 5 H;r anion. 
Still, the validity of this rule is rather limited, even if 
only monocyclic compounds are considered. The point is 
that far from all cyclic polyenes can have a planar non- 
stressed configuration. Whereas the rings of small size 
with (4w + 2) n-electrons are aromatic (the lengths of 
the chemical bonds in such systems being equal and the 
derealization energy contributing appreciably to the 
formation enthalpy), the systems with larger rings pos- 
sessing (4n + 2) n-electrons are no longer planar and have 
no aromatic properties. 
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However, when carbon atoms number 14, 18, and more, 
a molecule becomes a planar aromatic system again. 
Many years of persistent research had been spent before 
the compounds C 14 H 14 , C 18 H 18 and others were obtained. 
The authors called them [14] and'[18]annulenes respective- 
ly. [18]annulene is the first (after benzene) cyclic poly- 
ene with (An + 2) K-electrons, whose planar conforma- 
tion is stable at room temperature. It is less stable than 
benzene, but is more stable than the corresponding acy- 
clic polyene (nonaene): 

CH 2 =CH-(— CH=CH— ) 7 — CH=CH 2 . 

As to the [JVlannulenes with n < 18, they are not 
planar due to the repulsion of hydrogen atoms located 
inside the ring. In the case of [14]annulene the planar 
conformation is stable only at t = — 60°C. 

However, hydrogen atoms are not present inside the 
cycles such as the one shown below: 

* > 

In this and similar compounds the acetylene bond is 
supposed to donate only two n -electrons to the conjugated 
system while the other n-bond is located in the plane of 
the molecule and does not participate in the conjugation. 
Consequently, this compound satisfies the Hiickel rule for 
n = 4. It indeed possesses aromatic properties. 

Anti-aromaticism. When a cyclic polyene system is 
studied it is important to know whether this system is 
nonaromatic, i.e. not stabilized by conjugation and suf- 
ficiently reactive due to the internal tension and other 
causes, or destabilized by conjugation, i.e. the cyclic 
derealization increases the total energy of the system. 
In the latter case the molecule is called anti-aromatic. 

Here are typical examples of anti-aromatic systems: 
cyclobutadiene, a cyclopropenyl anion, a cyclopentadi- 
enyl cation, and others. 

As new cyclic conjugated systems were discovered the 
concept of aromaticity got more complicated and the 
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simple Hiickel rule became insufficient. Some liew, more 
general criteria of aromaticity were obviously necessary. 
Some of these criteria will be discussed below. 



In Search of Aromaticity Criterion 

Dewar's theory. In recent years more and more popular 
in organic chemistry becomes the theory of perturbation 
molecular orbitals (PMO) worked out by the American 
chemist M.J.S. Dewar. Using that theory, many com- 
plex problems of organic chemistry can be tackled, includ- 
ing the problem of aromaticity. Since it is impossible 
to describe this theory here, we are citing only some of the 
conclusions concerning the cyclic polyenes, drawn from it. 

Let us consider the interaction of two odd alternant 
hydrocarbons R and S: 



Such systems possess nonbonding MOs (Fig. 19) which 
on the formation of the molecule RS produce two MOs one 



of which, the lower one, is filled with two electrons. As was 
shown by Dewar, in the first approximation the change 
in the energy of the system (R + S) on the formation 
of RS is only due to the interaction of nonbonding MOs. 



R + S-+RS. 




R 



Fig. 19 
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If R and S are joined by means of the atom r in R and 
the atom s in S 

©~© - (JHD 

the indicated change of one-electron energy is determined 
by the following formula: 

where a r0 and b s0 are the AO coefficients of the atoms r 
and s in the expression for a NBMO. 

Then the change in the total n-electron energy will be 
twice as great: 

A£„-2a r0 6 s0 p rs , (29) 

since two electrons lost some of their orbital energy. 

If the two odd ASs react at two "points" of molecule 
each 



am 



the change in the total n -electron energy (in this type 
of bonding) is equal to AE„ for every two interacting 
points of the initial molecules: 

A£ n = 2a r0 & s0 p rs + 2a (0 & uo p,, t . (30) 

Let us consider a benzene molecule as an example. We 
may picture benzene resulting from bonding two allyl 
radicals: 

H 2 C^~^CH 2 
CH 

According to Eq. (26) the NBMO coefficients of an 
allyl radical are equal to 

a^^a a = i/V2. 
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Now we can easily obtain the value of which in 
the case of benzene equals: 




A£„ = 2a 2 P + 2a 2 p = 2p. 



For a molecule of 1,3,5-hexatriene: 




A£„ = 2a 2 p = p. 



Whence it is clear that the formation of benzene is 
a more favourable process (in terms of energy) than the 
formation of hexatriene because in the former case the 
total energy of the system decreases by a greater value. In 
other words, benzene is more stable than the correspond- 
ing acyclic triene. 

Dewar suggested the general method of assessing rela- 
tive stability of even alternant systems. To do this, a mo- 
lecule of such a system is to be mentally broken up into 
the methyl radical CH3 and the "residue" which is an 
odd AS. 

Then the "methyl" coefficient in Eqs. (29) and (30) is 
equal to unity, and the coefficients of the "residue" can be 
found by the method described on p. 59. Eventually we 
come to the same conclusions as before: benzene is more 
stable than hexatriene (Fig. 20a) and consequently is 
aromatic while butadiene is more stable than cyclobutadi- 
ene (Fig. 20b) and consequently the latter is anti-aroma- 
tic. 

Note the signs of the coefficients. We can conclude from 
what was said above that the higher or lower stability of 
a cyclic polyene as compared to an acyclic one depends on 
the combination of signs of the coefficients at the ends of 
the "demethylized" compound. If the signs are identical, 
the even AS is aromatic due to cyclic stabilization; if the 
signs are different, the system is anti-aromatic due to 
cyclic destabilization. Hence, the Hiickel aromaticity 



72 



rule can«be easily formulated: an odd AS with identical 
signs of NBMO coefficients contains (in + 1) n-centres 



(a) 



o -A 

CH. 



c 



AE„=2b-t-J3=2bfl 
A E„ =2bfi + 2bfi =4bfi 



a 

~ d CH, 



AE„=2afl 

AE„ =2o/J + 2(-a)p=0 



Fig. 20 



plus one ax-centre at the methyl radical, i.e. the total 
number of ji-electrons must be (in + 2). 

Polycyclic systems. Dewar has extended his method to 
include the polycyclic polyenes. Fig. 21 illustrates the 



decapentaenes 

&E„=2dfi 



(-*) 



1-phenylbutadiene \yo[-annulene 
AE„ = 4<tp 



or 




azulene 



naphthalene 



Fig. 21 



binding of methyl with nonatetraenyl resulting in the 
formation of various systems, both closed and with open 
chain. The highest stabilization occurs in the formation 
of naphthalene which actually possesses aromatic proper 
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ties, although not to the same extent as benzene. Similar- 
ly, it is easy to show that pentalene and heptalene are 
anti-aromatic (Fig. 22). 




AE„=2fl(a-a+o)=0 




AE„=2J3(a~a + 0)=0 

Fig. 22 

Aromaticity according to Dewar. The general Dewar 
rule can be formulated as follows: 

To determine the character of an even alternant 
cyclopolyene, one has to consider an odd alternant 
system obtained from the initial cycle by deleting 
one carbon atom. The binding of this odd AS to 
the methyl group results either in an even cyclic 
AS or in an even AS with one of its rings open. 
A cyclic system is more stable than any of the com- 
pounds with open chain, i.e. is aromatic, only when 
the NBMO coefficients of the respective atoms in the 
odd system have identical signs. Otherwise, the 
cycle is anti-aromatic. 
The Dewar aromaticity rule just cited holds if a consid- 
ered cycle possesses (4n + 2) n-centres and does not if 
the number of n-centres equals 4n, whence follows ... 



74 



...DewA- Hiickel aromaticity rule. In a polycyclic al- 
ternant hydrocarbon a cycle with (An + 2) n-centres 
is aromatic whereas a cycle with An jt-centres is 
anti-aromatic. 

Dewar also formulated some other rules according to 
which various cyclopolyenes can be classified with aro- 
matic compounds. However, we shall not dwell on them 
here, addressing the readers to the excellent monograph 
by M. J. S. Dewar and R. C. Dougherty*. 

Apparently, it became clear from the previous chapters 
that the Dewar theory is essentially topological although 
it does not use explicitly the terminology and methods of 
the graph theory. 

Modulo 4 rule. Thus, the Hiickel (An + 2) rule predicts 
that (An) annulenes are thermodynamically much less sta- 
ble than (An + 2) annulenes (e.g. compare benzene and 
cyclobutadiene). I. Gutman and N. Trinajstic advised to 
distinguish this rule from another regularity appearing 
also in the case of polycyclic conjugated systems, namely 
that the cycles containing (An + 2) ji-electrons are stable 
while those with An n-electrons are not (e.g. compare cyclo- 
pentadienyl anion and cation). 

Every cyclic conjugated structure can be characterized 
by its topology in terms of the ring mode. The mode of 
a ring can be expressed by one of the numbers from to 3. 
Thus, a ring system of mode corresponds to a An cycle 
(An = modulo 4), rings of mode 1 to (An + 1) systems 
(An + 1 = 1 modulo 4), etc. Since the Hiickel rule can be 
generalized to include polycyclic conjugated systems, 
rings of mode 2, i.e. (An + 2) systems, produce the stabi- 
lizing effect whereas rings of mode produce the destabi- 
lizing effect (cf. cyclobutadiene, benzocyclobutadiene, 
and biphenylene). The size of a ring is also important. As n 
becomes sufficiently large, the difference between the 
rings of mode 2 and mode gradually disappears. 

Although rings of modes 1 and 3 do not stabilize or 
destabilize a molecule, their presence in the polycyclic 
conjugated systems manifests itself in their red-ox behav- 
iour. Hence, the systems containing rings of mode 1 can 
be easily transformed into the anionic form while mole- 



* M.J.S. Dewar, R. C. Dougherty. The PMO Theory of Organic 
Chemistry. A Plenum/Rosett Ed., N. Y., 1975. 
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cules containing a ring of mode 3 can be easily "oxidized, 
that is transformed into the cationic form. 

Topological Ji-resonance energy. The concept of the 
ji-resonance energy, or, in other words, the derealization 
energy (see p. 24), of conjugated molecules has served for 
many years as a simple criterion of aromaticity. The 
it-resonance energy (RE) in {5 units is defined as 

A£n E = £ n (conjugated molecule) — 

— E„ (hypothetical reference structure) 

(see p. 18). 

The main difficulty with the RE concept is the hypo- 
thetical nature of the reference structure, its choice being 
somewhat arbitrary. There are many ways of defining RE 
(M. Dewar, C. F. Wilcox, and others). Here we shall brief- 
ly examine the topological resonance energy (TRE). 
TRE is a nonparametric index which is directly related 
to the topology of a molecule and is of great practical 
value in predicting the aromatic stability of an arbitrary 
n -network. 

The TRE can be simply expressed as follows: 

N 



where N is the number of vertices (atoms) in a graph (con- 
jugated molecule) G, gj is the occupancy number of the 
/-th MO, X) (J — 1, 2, . . ., N) are the roots (eigenvalues) 
of the characteristic polynomial of the given graph (mole- 
cule) G, or, in other words, the individual orbital ener- 
gies, and Xj° are the roots of the acyclic polynomial corre- 
sponding to this graph and called a reference polynomial: 

N 

P ac (G; x) = T a« c ^-«, 

n=0 

< c = 5 (-i) C(S \ 

S£S«c 

where S£ c are the acyclic Sachs graphs having exactly n 
vertices. 

The TRE values thus obtained are used to analyse and 
classify arbitrary conjugated systems: the molecules with 
TRE (per electron) > 10" 2 are considered to be aromatic; 
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if — 10'-* < TRE (per electron) < 10" 2 the molecules 
are nonaromatic; and when TRE (per electron) <—10~ 2 the 
molecules are classified as anti-aromatic (see below). 

© 

□ 

0.046 0.039 -0.003 -0.307 

Still another (obviously topological) approach to the 
problem of cyclic conjugated systems are proposed in 1971 
by the American scientists M. J. Goldstein and R. Hoff- 
mann. The next section will be devoted to it. 

Pseudosymmetry, Topology, Aromaticity 

How to make ends meet. Let us assume that each cyclic 
conjugated hydrocarbon consists of several polyene rib- 
bons connected in a certain way. Our objective is to find 
out how the way in which such ribbons are interconnected 
affects the stability of an organic molecule or ion com- 
posed of them. 

For one ribbon there is the unique mode of connecting 
its ends: 

, for example 

When a system comprises two polyene ribbons, the 
latter can be connected in two ways which are topological- 
ly nonequivalent: 




In the first case each end of one ribbon is connected only 
to one end of the other ribbon. Such a topology of connect- 
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ing the ends is called pericyclic. In the second case 
each end of one ribbon is connected to both ends of the 
other ribbon (the spirocyclic topology). 

Topology of connecting three ribbons is even more di- 
versified. Here are the most important examples: the peri- 
cyclic connection of ribbons 




the longicyclic connection of ribbons ("keg") 




and the laticyclic connection of ribbons ("ladder") 




Permitted interactions of ribbons. To simplify the 
further discussion, let us impose some restrictions on the 
ribbon interaction mode. First, we assume that the rib- 
bons interact only with their ends. Second, the ribbon 
twisting angle should not exceed 90°. 
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And'ijnally, the ends of each ribbon interact quite 
identically, i.e. with an equal number of other ends and 
via formation of bonds of the same type, either a or ji. 

Pseudosymmetry of ribbon's orbitals. In our following 
discussion we shall be looking for symmetry of ribbon's 
orbitals in the case of reflection in the plane *P (Fig. 23). 



>iiS VA 
Fig. 23 

Evidently, this plane is not a true symmetry plane be- 
cause the specific geometrical properties of a ribbon are 



<PA 



GO 



n l 



2 3 4 

Fig. 24 
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disregarded. In particular, the ribbon twisting is ne- 
glected as well as the presence of various substituents at 
its ends. Here the term "reflection in a plane" do es no t 



imply a mirror reflection of an object. Actually we deal 
here with a more general problem of examining whether 
an object gets inverted on transition "beyond the mirror" 
or not. In other words, the question is whether the signs 
of terminal />„-AOs of the ribbon change on reflection in 
the plane V or not. Depending on that, we deal with 
pseud osymmetric (WS) or pseudoantisymmetric^-^MOs 
respectively. The plane W is called a pseudosymmetry plane. 

Alternating pseudosymmetries. Looking at Fig. 24 show- 
ing the relative energy positions of MOs for the polyene 
series, one can easily notice that the pseudosymmetries 
of these MOs alternate, the lowest energy orbital being 
of the WS type. 

Classification of ribbons. Now let us pass from the 
classification of molecular orbitals of polyene fragments 
to classifying the ribbons. It can be done in different ways. 
Sometimes it is sufficient to indicate only the number 
of it-electron centres (n) and the charge (z) of a polyene, 
denoting each ribbon by the symbol n z . Then, for exam- 
ple, the fragment 



will be denoted by the symbol 4°, and the fragments 



© and 




by the symbols 3" and 1 + respectively. 

But the ribbons can be (and, for our purposes, must be) 
classified differently, taking into account the population 
density and pseudosymmetries of the highest occupied 
MOs (HOMO) and the lowest unoccupied MOs (LUMO). 
Then all polyene ribbons get subdivided into four groups 
(see the upper part of Fig. 25). To avoid the tiresome 
drawing of the MO layouts given above, we encode each 
layout by a certain number. Goldstein and Hoffmann ad- 
vised to do it as follows: to each ribbon we ascribe the num- 
ber (n — z) mod 4, where the arithmetic operation mod 4 
stands for subtraction from (n — z) of the maximum 
multiple of four, i.e. Aq (q = 0, 1, 2, . . .)*. In other words, 

* Note that (n — z) is equal to the number of Ji-electrons in a 
system. 
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(« — z) mod 4 is the remainder of division of (n — z) 
by 4. Obviously, Aq must be less than or equal to (n — z). 
Therefore, if (« — z) < 4, the subtrahend must be zero 
(q = 0). Thus, for example, the butadiene fragment (4°) 



lumo 


VA — <PA — 


— 


IPS — 


HOMO 








(n-z) mod 4 : 


1 4 


3 





Examples : 


2 + 2° 


2" 


2" 




3 ++ ,3~~ 3 + 


3° 


3" 




4" 4 ++ ,4" 


4t 


4° 




5° 5" 


5 ++ ,5" 


5 + 




6 + 6° 


6" 






7 ++ ,7" 7 + 


7° 






Fig. 25 







has a corresponding number (4 — 0) mod 4=4 (mod 4) = 
= 4 — 4 = and the butadienyl anion (4~) is associated 
with the number (4— (—1)) mod 4 = 5 (mod 4) = 1, 
and so on. 

It is easy to surmise that the expression (n — z) mod 4 
takes only the values 0, 1, 2, 3, each of which can be 
brought into a one-to-one correspondence with one of the 
layouts given earlier (Fig. 25). 

Acyclic connection of ribbons. Thus, we have finished 
the preliminary work and can set to connecting ribbons. 

First, we shall connect only two ribbons and only at 
one location. The result will be a new acyclic polyene 
fragment. Let us call it intermediate. If the interacting 
orbitals of the initial ribbons had an identical pseudo- 
symmetry, e.g. ¥5, they would combine to form two 
MOs of the intermediate ribbon with different pseudo- 
symmetries, WS and WA. In this case, the energy 
of the intermediate MO is lower, and the energy *¥A is 
higher, than the energies of the initial WS MOs (Fig. 26a). 
The same is true when two MOs interact (Fig. 266). 



6-01102 
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But if the interacting MOs have different pseudosym- 
metries, the energy WA of the intermediate MO will be 
lower than the initial orbital energies (Fig. 26c). Note 



(a) 

>i>s + </>s=:^ 

(b) 



<PA 



(c) 



<M4 + <M4=< ^ A 

, - 

VA 

Fig. 26 

that the interacting terminal orbitals of all intermediate 
MOs possessing higher energies are in antiphase and there- 
fore do not overlap. That is why they are located higher 
than other MOs on the energy scale. For low-lying MOs 
the opposite is true. 

A ring is completed. Let us close the cycle. In the 
process, the low-lying intermediate MOs must also get 
"closed", that is, their terminal p electron clouds must 
get overlapped. However, this is possible only when the 
MOs of the initial ribbons are of identical pseud osymmetry 
(Fig. 26). 

Besides, on closing the cycle the energy of intermediate 
MOs changes in accordance with ,their pseudosymmetry: 
the energy of WA MO increases and that of WS MO de- 
creases (Fig. 27). 

We avoid using here the terms "bonding", "antibond- 
ing", and "nonbonding" MO because in the framework of 
the Hiickel method the HOMO of a ribbon of a certain 
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(e.g. 0) type may be bonding (for systems 4°, 8°), nonbond 
ing (3", 7~), and even antibonding (2 , 6 ). Surely, the 
same is also true for a LUMO. 

Stability of pericycles. Hereafter to simplify our pre- 
sentation we shall confine ourselves to the cases where 
only the polyene ribbons of even types, and 2, interact. 



Ribbon type . 
<PS 



Fig. 27 



tyA 



Accordingly, when examining the pericycles formed by 
two ribbons, it is sufficient to consider the following 
three cases of ribbon interaction: + 0, + 2, and 2 + 2. 

Suppose the interaction is between two ribbons of 
type 0. Then as is seen from Fig. 28a, the pericycle forma- 
tion lacks effective stabilization: two electrons increase 



(0.0) (0,2) 2 (2,2) 2 




(a) (b) (c) 

Fig. 28 



their energy while two others decrease. An example of 
such a nonstabilized system is cis-9,10-dihydronaphtha- 
lene (4° + 4°) which was first obtained in 1963: 




The (2 + 2) interaction (Fig. 28c) is similar to that 
just described. As an example, we can indicate the mole- 
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cule of norbornadiene (bicyclo[2,2,l]-hepta«2,5-diene) 
relating to type (2° + 2°) and being a fairly reactive 
compound: 




Another example is Dewar's benzene (bicyclo[2,2,0]- 
hexa-2,5-diene) obtained in 1963: 



The latter is an unstable compound. When kept at room 
temperature it gradually turns into common benzene, 
with only about half of the initial quantity left two days 
later. 

And only the pericyclic interaction of two ribbons of 
different types, i.e. the (0 + 2) interaction, brings about 
a stabilized system (Fig. 286). Once again we draw the 
reader's attention to the fact that the relative positions 
of MOs in Fig. 28 are shown as if only the interaction be- 
tween MOs of identical pseudosymmetry existed, in 
accordance with what was stated earlier. Henceforth we 
shall be denoting the type of the obtained cyclic hydro- 
carbon by the symbol (k, I), where k and I are the types 
of the initial forming polyene ribbons. 

An example of a relatively stable pericycle of type 
(0, 2) is bicyclo[4,2,l]-nona-2,4,7-triene 



whose derivative 9-methylbicyclo[4,2,l]-nona-2,4,7- 
triene-9-ol was first obtained in 1963 from aromatic 
cyclooctatetraene dianion. 

Bicyclo[4,2,l]-nona-2,4,7-triene is quite stable ther- 
mally. It sustains heating up to 450°C and only under 
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ultraviolet radiation undergoes various regroupings of 
atoms, forming primarily tricyclononadienes. 

Some compounds of types (2, 2) and (0, 0) regroup com- 
paratively easily into (0, 2) pericycles. An example is 
the spontaneous transformation of unstable 1,6-methane- 
cyclodecapentaene into double norcaradiene: 




The pericycle stability results obtained are summarized 
in Table 1. 

Spirocycles. Now let us turn to the spirocyclic topology 
of connection of two ribbons. As in the previous case, we 
shall take into account only the interactions of orbitals 
possessing identical pseudosymmetry. The new topology, 
however, calls for imposing new restrictions on the ribbon 



Table 1 



Num- 
ber of 
poly- 
ene rib- 
bons 


Cycle 
topology 


Stabilized 
cycle 
type 


Number of 
n-elec- 
ons in 
reach 
ithon 


Total 
number 
of n-elec- 
trons 


Example 


2 


Pericyc- 
lic 


(0, 2) 


Aq 
+ 2 


47 + 2 


Bicyclo[3,2,lJ- 
octadienyl 












anion 














<? 


7 












3" + 2° 




2 


Spiro- 
cycles 


(0, 2) 


Aq 
4? + 2 


4 ? + 2 


Spiro[2,4J- 
hepta-1,4,6- 
triene 

CK 

4° + 2° 





interaction modes. The interaction of two MPs of 
type involves neither essential lowering nor increasing 
of their energies. Therefore, when studying stable spiro- 
cycles special attention should be paid to WA MOs. 

Again, the three cases of ribbon interactions must be 
considered: + 0, + 2, and 2 + 2. Similar reasoning 
enables us to come to the following conclusion: only spi- 
rocycles of type (0, 2) are stable. 

Unfortunately, experimental data on spirocyclic sys- 
tems are insufficient so far and the number of examples 
is thus restricted (Table 1). 

Longicycles. As was noted before, the longicyclic 
topology is inherent in compounds possessing at least 



three polyene ribbons. Let us derive the conditions for 
a stabilization of longicycle. To do this, we shall consid- 
er Fig. 29 illustrating various cases of terminal-orbital 
overlappings of three polyene ribbons constituting a lon- 
gicycle. Let us "forget" for a white about the third (up- 
per) ribbon. Then the remaining pericycle can be charac- 
terized by indicating the pseudosymmetry of terminal 
orbitals of constituent ribbons with respect to the planes 
W 1 and W 2 . If the orbitals of the two pericyclic ribbons are 





Fig. 29 
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of the type WAA or WSA , the ends of the third ribbon can 
effectively interact with these orbitals and all three rib- 
bons become thus connected. In this case the third ribbon 
must be pseudosymmetric relative to the plane ^ if the 



o 0+(o,o) (o.o) (0,0) + 2 2 
fy-<t>SA ■ 



<PSS —h— 

M — +t-</MS 4f — H- 




+ (0,2) {0,2) (0,2)+2 




±<PAS *. 




IPS 



0+(2,2) (2,2) (2.2)+2 
-<PAA 




0A 



'PA *iXV \)*<t>S 



-VAS 



<PSA + 



Fig. 30 



two other ribbons are of the type WSA , and pseudoanti- 
symmetric relative to the same plane if the pericycle is 
of the type ¥.4.4. In other words, the effective longicyclic 
connection of three polyene ribbons results only in the 
case of the (WA + WAA) and (VS + VSA) interactions 
(see Fig. 29). 
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As before, we shall consider only the interaction of 
ribbons of the even type, i.e. we confine ourselves to 
the cases 

+ (0 + 0) 
+ (0 + 2) 
+ (2 + 2) 

and 

(0 + 0) + 2 
(0 + 2) + 2 
(2 + 2) + 2 

The parentheses signify that we consider first the 
orbitals of two pericyclic ribbons and then the change 
of orbital energies after their interaction with the third 
ribbon is taken into account. The corresponding molecu- 
lar orbital diagrams are shown in Fig. 30. 

Studying these diagrams, one can come to the follow- 
ing conclusion: stabilization is natural for the longicycles 
of types (0, 2, 2) and (0, 0, 2)*. 

The longicycles of other types, as, for example, bar- 
relene (type (2, 2, 2)) 




and 7-norbornadienyl anion (type (2, 2, 2)) 




are nonstabilized. 



* Unquestionably, this conclusion is independent of the order 
according to which we "switch on" the ribbon interactions: (0+ 2) + 
+ 2 or + (2 + 2), etc. 
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(0, 2, 2) 



At the, same time, 7-norbornadienyl cation of type 




is quite stable. In 1968 data were obtained on the existen- 
ce of a stable anion 



Laticycles. The studies of laticycles carried out by 
Goldstein and Hoffmann disclosed that compounds of the 
types (0, 0, 2), (2, 2, 0), (0, 2, 0), and (2, 0, 2) must be 
stable, the two last types being especially so. Unfortunate- 
ly, there are no sufficiently reliable experimental data 
so far to confirm their conclusion. 

Tips on terminology. Earlier, when characterizing one 
or another system, we used the adjectives "stabilized" 
and "nonstabilized". 

It was noted in Chapter 3 that the stability of a chem- 
ical compound is a very complex concept incorporating 
many different aspects. Therefore, in each specific case it 
is necessary to explain what is meant. Here the expression 
"stabilized cycle" signifies that the formation of such 
a cycle from polyene ribbons results in lowering of the 
energy of corresponding n-MOs. Otherwise, the cycle is 
said to be destabilized. In other words, we compare the 
stability of a cyclic hydrocarbon to that of isolated 
polyene ribbons. In doing this, we look only into the 
n-electron energy of the system. 

What was said above clearly indicates that the destabi- 
lization of the cycle must lead, as a rule, to higher re- 
activity of the compound. The opposite statement, how- 
ever, is not generally true since even stabilized systems (in 
the sense indicated) can be rather reactive either due 




0- 



to the presence of electric charges or due to Unpaired 
electrons in them, or owing to some other causes. 

New definitions. Thus, the stabilization of cyclic 
conjugated systems depends to a great extent on their 
topology. In fact, the hydrocarbons formed by three 
polyene ribbons and related to the type (2, 2, 2), i.e. pos- 
sessing Aq + 2 n-electrons, are stabilized in the case of peri- 
cyclic topology* (e.g. benzene) whereas the correspond- 
ing longicycles and laticycles are destabilized. However, 
in the case of "three-ribbon" hydrocarbon of the type 
(2, 2, 0) the oppositejstatement is true. Hence, it becomes 
possible to predict certain rearrangements proceeding 
along with the transformation of the system's topology, 
for example: 




pericycle Lcngicycle 



Goldstein and Hoffmann suggested the new terms: bi- 
cycloaromaticity and antibicycloaromaticity. The bicy- 
cloaromatic compounds are those cyclic conjugated hyd- 
rocarbons which contain 4g jt-electrons and stabilized 
in the shape of laticycles and longicycles. The cyclic 
conjugated systems which possess 4g + 2 n-electrons 
and are destabilized in the shapes indicated were ac- 
cordingly proposed to be called antibicycloaromatic. 

Future will show how useful is such a classification, 
but it is indisputable that the cycle topology proves to 
affect essentially the properties of the cycle. 



* Although earlier we did not consider the conditions lor sta- 
bilization of pericycles formed from three ribbons, the result in- 
dicated can he easily derived by the reader himself. 



Chapter* Five 



Molecules 
Without 

Chemical Bonds 




So far we discussed the application of the graph theory 
in studies of electronic structure of molecules. In this 
chapter we shall describe some topological concepts that 
appeared in stereochemistry in connection with the 
investigations of a new remarkable type of chemical 
compounds, catenanes, rotaxanes, knots. 



Topological Isomerism and Topological Bonding 

Let us take two wire rings which are not connected to- 
gether. Nothing prevents us from arranging them arbitrar- 
ily relative to each other. Let us pass one of the rings 
through the other. Now we cannot arrange the rings arbi- 
trarily any more. 

A similar situation is observed in cyclic molecules. 
"Passing" one of them through another, we can obtain 
a new configuration in which two cycles, though nonbond- 
ed chemically, constitute a single whole. Such a bonding 
is called topological, and the composite chain molecule 
is called a topological isomer of the sum of two individual 
constituent molecules. The compounds of this kind are 
classified with catenanes (from Latin catena, chain) 



Let us consider another way of connecting molecules 
without chemical bond formation. If through a macrocy- 
clic molecule passes a linear one with bulky terminal 
groups, the system thus obtained represents a topological 
isomer of the type "dumbbell in a ring" (Fig. 316). Such 



(Fig. 31a). 
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compounds are called rotaxanes (from Latin vords rota, 
a wheel, and axis). 

Still another type of topological isomerism is "molecule 
in a cage". Formally, here belong the inclusion compounds 




Fig. 31 



(clathrates), carbide organ ometallic clusters (cryptates) 
(e.g. an alkali metal ion inside a macrocyclic ligand), and 
other systems. However, there is apparently a chemical 
interaction in them between a particle trapped inside the 
system and its surroundings. 

Also known are other types of topological isomers 
which draw the researchers' attention, but they exist only 
on paper so far. Catenanes and rotaxanes are the only 
topological structure classes which have been investigated 
experimentally. 

Links of a Single Chain 

Cycle twining. The idea of obtaining molecules consisting 
of two cycles passing through one another was discuss- 
ed by German chemist R. Willstatter as early as the 
beginning of this century. However, the first efforts to 
synthesize catenanes failed. To pass one ring through 
another, it was necessary to obtain a large monocycle 
composed of at least 20 atoms. Such macrocycles were 
obtained at the middle of the 50s. In 1956-58 the problem 
of catenane synthesis was actively attacked simultaneous- 
ly by several research teams: one headed by A. Liittring- 
haus at Freiburg University, another by F. Cramer at 
Heidelberg, and still another by R. Kohler and B. Diet- 
rich at Tiibingen.... No success again. 

Statistical synthesis. As far back as in 1952 an idea was 
ventured concerning an accidental formation of a catenane 
during the polymerization reaction, e.g. in polyphospho- 
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nitryl chloride. In I960 E. Wasserman proved reliably 
the catenane formation for the first time. He reasoned 
that the cyclization of linear molecules in the presence 
of cyclic ones provides a chance (although a very scarce 
chance) for the chain to pass through the ring, 



as well as for the mutual "interception" of linear mole- 
cules and the eventual formation of a catenane: 



The probability of free intertwining grows with the 
length of molecules, but such molecules are very bother- 
some to work with because their solubility decreases dras- 
tically as the chain gets longer. Wasserman subjected 
diether ROOC — (CH 2 ) 32 — COOR to acyloin {conden- 
sation in the presence of partially deuterated cyclotetra- 
triacontane (C 34 H 63 D 5 ). The latter (not deuterated, of 
course) was obtained back in 1926 by the famous Swiss 
chemist L. Ruzicka during the analysis of natural aro- 
matic musk substances. 

Here is the Wasserman synthesis chart: 



+ 






+ C 34 H 63 D S + 



+ acyloin 
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Next, the obtained mixture was separated in u chroma- 
tographic column using alumina and silica gel. In the 
process, C 34 li 63 D 5 was completely removed by washing 
the column with pentane, and acyloin was washed off by 
an ethermethanol solution. However, even after repeated 
washings the deuterated cyclic compound C 34 H 63 D 5 re- 
mained in the column. And only after the decomposition 
of the acyloin cycle 



was the Ruzicka hydrocarbon at last drawn off the column. 

These data indicated that the acyloin ,cycle is inter- 
twined with C 34 H 63 D V Initially Wasserman recovered 
5.66 mg of catenane oil, the total yield being 0.0001 per 
cent. Subsequently Wasserman and coworkers developed 
other methods of statistical synthesis of catenanes. 

Directional synthesis. In April of 1964 in the laboratory 
of the Freiburg University G. Schill and A. Liittringhaus 
succeeded in the directional synthesis of a catenane ac- 
cording to the following chart: 

ci /(CH a ) 12x 




+ HOOC — (CH 2 ) 3J — COOH 





HO 



✓(CH,),^ 
(CHs) a5 C=0 



5,0 *~ 
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Being rather complicated (twenty stages!*), this proce- 
dure nevertheless is reliable, ensuring the correct identi- 
fication of a catenane, and provides a fair 20 per cent 
yield. 

In 1969 G. Schill and C. Ziircher synthesized a catenane 
consisting of three rings passing through one another: 




Ac 



In addition to catenanes, the rotaxanes were also ob- 
tained through both statistical and directional syntheses. 
One of them is shown below: 




Characteristically, R have to be bulky groups. 

Catenanes in nature. In 1967 the American biochemist 
G. Vinograd found natural catenanes. Using an electron 
microscope, he obtained photographs showing clearly 
intertwined catenanes up to a [7]-catenane. The size of 
individual macrocycles reached 5 um. Later it was 
demonstrated that catenane DNAs are common in nature. 

Syntheses of tomorrow. The search for various topo- 
logical isomers in nature and the efforts of their laboratory 
synthesis continue. Many hypotheses and reaction charts 



Naturally, only few stages are indicated here. 



were suggested, and the properties of expected campounds 
predicted. 

For many years chemists try to obtain twisted molecules 
resembling the Mobius band. The feasibility of such 
systems was discussed for the first time by H. Frisch 
and E. Wasserman as early as in 1961. 

In some cycloolefin transformation reactions in the 
presence of transition metal compounds certain intra- 



molecular arrangements can be expected with the forma- 
tion of knotted and twisted molecules (necktinodanes and 
plectanes) (Fig. 32). The synthesis of such compounds 
is an accomplishment of tomorrow. 




IKXXX 



Fig. 32 



chapter six Topological 

Concepts 
in Inorganic 
Chemistry 



Until now we considered the topology of various or- 
ganic molecules. In inorganic chemistry the topological 
methods are not so popular. Therefore, in this chapter 
we shall confine ourselves to only two topics, the topolo- 
gical theory of borane structure and the discussion of the 
interrelationship between topology, symmetry, and ele- 
ctronic structure of coordination compounds. 



Boranes: Structure, Bonding, Topology 

Chemistry of boron hydrides (another name, boranes) is 
now one of the most challenging and rapidly growing area 
of inorganic chemistry. Boranes are employed in many 
fields of man's activity, primarily as missile fuels and 
in cancer research. 

In addition, boranes are of great interest for the chem- 
ical bonding theory on which we shall dwell further. 

First achievements. The founder of borane chemistry is 
rightfully held to be the German chemist A. Stock al- 
though the first boron and hydrogen compounds were ob- 
tained back in 1881. But, according to N. V. Sidgwick, 
"All statements about the hydrides of boron earlier than 
1912, when A. Stock began to work upon them, are untrue". 

During 25 years of persistent work Stock and his col- 
leagues managed to obtain (and determine the physical 
and chemical properties of) the following boron hydrides: 
B 2 H 6 , B 4 H 10 , B 5 H 9 , B 5 H U , B 6 H 10 , and B 10 H 14 . With the 
exception of diborane, these compounds were obtained 
through the interaction of magnesium boride Mg 3 B 2 with 
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diluted acid. Diborane B 2 H 6 was obtained by thermal de- 
composition of higher boranes. These days, as a rule, 
a different approach is used: first, diborane is obtained by 
the reaction 

In ether 

3NaBH 4 + 4BF 3 > 3NaBF 4 + 2B 2 H 6 

or by some other method, and then the higher boranes are 
produced from diborane. 

Deceptive resemblance. For a long time chemists were 
occupied with the problem of the structure of diborane and 
other boron hydrides. In particular, the similarity in 
empirical formulas of ethane and diborane: C 2 H 6 and B 2 H 6 
attracted attention. That is why many chemists initially 
accepted that diborane has an ethane-like structure. Ma- 
ny..., but not all. The first who had his doubts was the Ger- 
man chemist W. Dilthey. In 1921 at the meeting of the 
Erlangen Chemical Society he presented his paper "On 
Structure of Water" in which he suggested the following 
structure of dimeric water molecules (H 2 0) 2 : 



/ \ / 



At the same time Dilthey noted that B 2 H 6 can also be 
described by a similar model: 

n b' ;b x 



H 

However, Dilthey's hypothesis was soon forgotten, and 
the idea of similarity with ethane prevailed. Only in the 
40s the arguments concerning the diborane structure re- 
newed. 

As early as 1940 the American chemist F. Stitt demon- 
strated that in B 2 H e the barrier hindering the internal ro- 
tation about the hypothetical B — B bond is much higher 
as compared to ethane. Then some more experimental data 
appeared to promote the hydrogen bridge structure. Dil- 
they's model was again taken up in the publications by 
B. V. Nekrasov, Ya. K. Syrkin, M. E. Dyatkina in the 
USSR, and some foreign chemists. 
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The final proof of the presence of hydrogen bridges in 
boranes was provided in 1948 for decaborane B 10 H 14 and 
in 1951 for B 5 H 9 and B 2 H 6 . 

The theoretical description of similar structures in 
terms of a three-centre bonding conception (see below) 
was first suggested in 1943-47 by H. C. Longuet-Higgins 
and R. S. Mulliken. A great contribution to the investi- 
gation of structure and properties of boranes was provided 
by the American chemist W. N. Lipscomb who was 
awarded the Nobel Prize for his efforts in 1976. 

Deficiency without deficiency. Let us turn now to the 
electronic structure of boranes. Usually they are classified 
with electron-deficient compounds because the number of 
possible two-centre interactions exceeds the number of 
valence-electron pairs. 

In a molecule B 2 H 6 each boron atom forms two common 
(two-centre two-electron) covalent bonds with terminal 
hydrogen atoms. As to the bridge bonds B — H — B they are 
considered as triangular fragments in which the formation 
of three MOs involves one valence orbital of each atom 
(Fig. 33a). Each fragment contains two electrons occupy- 
ing the bonding MO W v That is how a very strong two- 
electron three-centre bond BHB appears which is strong- 
er than the conventional bond B — H. Thus, there is 
actually no deficiency of electrons. That deficiency exist- 
ed only in the earlier conceptions of chemical bonding 
which was always pictured as a result of two-electron two- 
centre interactions. 

In higher boranes, in addition to the bridges BHB, 
there are also two types of fragments with three-centre 
bonds: the "open" bridge B and the "closed", or 



"central", structure B. 



B B 

The molecular orbital charts corresponding to these 
fragments are illustrated in Fig. 336 and c. 

Lipscomb's topological theorem. In 1954 W. Lipscomb 
with coworkers suggested a so-called topological theory 
of boranes which he continues to improve even now. 
On the basis of that theory it became possible to predict, 
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and then tfo synthesize, some boranes and their derivatives 
(B 7 H 15 , B 9 Hi5, B U H 17 etc.). Lipscomb's topological the- 
ory is based on the balance equations derived by him and 
relating the borane composition with the number of cer- 
tain structural elements in borane. Lipscomb distin- 
guishes the following basic structural units: 



B 



B— H B B B— B B B and B ,x^ B 

As an example, let us consider neutral boranes with 
composition B p H p+g . Suppose the number of hydrogen 
bridge atoms is denoted by s, the number of terminal BH 
or BH 2 groups by x, the number of three-centre bonds BBB 
(both open and central) by t, and the number of two- 
centre bonds B— B by y. 

Then for hydrogen atoms the balance equation takes 
the form*: 

s + x = q. (31) 
Since each boron atom has four orbitals accommodating 
three electrons, the total number of three-centre bonds 
(BHB and BBB) is equal to the number of boron atoms in 
the molecule: 

s + t = p. (32) 

An example is a decaborane-(14) molecule whose topo- 
logical graph has the form: 




In this case s = 4, t = 6, and p = 10, and Eq. (32) is 
thus satisfied. 
And finally, the last balance equation 

t + y + q/2 = p (33) 

means that each boron atom in the structures in question 
is connected at least to one terminal hydrogen atom. As an 

* Make sure this equation holds, e.g. in the case of diborane. 
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example, the reader is advised to consider th-3 graph of 
decaborane-(14) shown above (t = 6, y = 2). 

For borane ions with the charge c (B p Hp +g+c ) the 
balance equations take the form: 

s + x = q + c, 

s + t = p + c, (34) 
t + y = p — c—±. 

The application of these equations can be illustrated 
for the case of the dianion B 10 H*;. According to our nota- 
tion, it should be written as B 10 H*~ +6 _ 2 . This ion is repre- 
sented by the following topological graph 

H 2 B B j B — BH.. 
H 

for which s = 2, t = 6, y = 3, x = 2, c = — 2, p = 
= 10, q = 6. 

From the Lipscomb theory one can get the topological 
graphs for any borane of assigned composition and thus 
predict new structures. Certain limitations should be 
imposed, though, on the quantities appearing in the bal- 
ance equations (31)-(34), for example, 

<?/2 <; p (or q). 

Besides, one can take into account the limitations due 
to the symmetry and metric characteristics of compounds 
(e. g. the limiting values of valence angles etc.)*- So, for 
borane B 10 H 14 (p = 10, q = 4) the following sets of 
the variables s, t, y, x (2sSC 4) are possible: 

styx 
4 6 2 
3 7 11' 
2 8 2 



* For example, the neighbouring boron atoms cannot be simul- 
taneously connected by central and bridge three-centre bonds BBB 
since in that case the valence angles would become very small. 
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Thus, tfye topology corresponding to the set of indices^ 
styx = 4620 is not the only possibility for a borane '•afcj 
composition B 10 H 14 . The problem of the existence of 
a specific structure can be solved both experimentally 
and theoretically, by the MO method. 

Weakness of the theory. Lipscomb's theory based on the 
assumption of localized three-centre bonds becomes inade- 
quate as the number of hydrogen and boron atoms in 
molecules grows, increasing dramatically the number of 
possible resonance structures that should be taken into 
account when composing the total wave function of the 
molecule. Inevitably, the number of possible topological 
types of the molecule also grows with the increase in 
p and q. 

It does not imply, however, that the topological theory 
has lost its significance to become second to the MO meth- 
od. It would be more appropriate to say that the molecu- 
lar orbital and topological approaches are mutually com- 
plementary. That it is really so can be seen from one of 
the recent papers by Lipscomb and his coworkers which 
is devoted to the study of unstable boranes B 2 H 4 , B 3 H 7 , 
B 3 H 9 , B 4 H 8 , and B 4 H 12 . 

New refinements. Using the MO method to calculate the 
optimal geometry of these compounds, it was established 
that the molecules B 2 H 4 , B 3 H 9 , and B 4 H 8 have vacant 
orbitals at a boron atom. To describe the unstable boranes 
with vacant orbitals, the known topological rules (the 
balance equations) had to be extended. The new equations 
for the systems BpH p+g take the form 

s + x = q, 

s + t = p — v, 

t+y = p—\, 

where v is the number of vacant orbitals. 

The transformation energy values of the boranes in ques- 
tion calculated by the MO method and by the semi-empi- 
rical equation proposed by the authors 

E = — A5v + 17 As + 28A* 
agreed to within 21 kJ/mol. 
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This result emphasizes once more that it pays to com- 
bine rationally simple methods of treating chemical com- 
pounds and complex quantum-mechanical calculations. 

Symmetry and Topology of Complexes 

Now let us turn to another aspect of graph-theoretical 
methods used in inorganic chemistry. We shall deal with 
the determination, even though approximate, of the se- 
quence of values of molecular-orbital energies in complex 
compounds. From the very beginning we should note 
that everything to be discussed can be attributed also to 
organic molecules possessing sufficiently high symmetry 
and comprising an adequate number of atoms at topo- 
logical^ equivalent positions. However, at present the 
approach described here is applied almost exclusively 
to high-symmetry inorganic compounds. 

Topology of interligand interactions. Consider, for 
example, a tetrahedral molecule AL 4 , where A is the cen- 
tral atom and L are identical a-ligands. The topological 
matrix of such a molecule can be written in the form 

(k 1 1 1 1\ 
1 I 
1 I, (35) 
1 I 
1 / 

where k is a certain parameter characterizing the differ- 
ence of the atom A from the ligands. Then the problem 
of determining orbital energies by the Hiickel method 
(since we treat our complex in terms of that method) 
reduces in the orthogonal basis (S = I) to finding the ei- 
genvalues of the topological matrix A. For small k we 
obtain: 

X i. 2, 3 = 0, 

However, some scientists believe that the matrix (35) 
has an essential drawback: topologically, it possesses re- 
dundant information, describing the topology of the mole- 
cule and, through the parameter k, the difference in the na- 



104 



ture of the central atom and the ligands. This drawback 
can be done away with in various ways. Thus, it was pro- 
posed to split the topological graph of the system AL 4 
into two subgraphs (Fig. 34) and then to examine the to- 
pological matrix of the subgraph G 2 . This procedure is 
equivalent to regarding only interligand interactions. 



L L 




G ==> G, © G, 
Fig. 34 

Accordingly, the order of increase of orbital energies ob- 
tained thereby is defined only by the ligand-ligand interac- 
tions. As to the interactions of the ligands with the cen- 
tral atom (exceeding in magnitude the interligand ones!) 
and the modifications in molecular-orbital charts, 
appearing when these interactions are taken into account, 
they can be calculated by the methods of quantum-me- 
chanical perturbation theory. Although this approach 
contradicts the adopted ideas of the relative intensity of 
ligand-ligand and ligand-central atom interactions, it 
turns out to be quite satisfactory for qualitative analysis 
of electronic structure of complexes. 

Sequence of orbitals. The topological matrix of the 
graph G 2 has the form: 

(0 1 1 1\ 
10 111 
110 1 J ' 
1110/ 

To find the eigenvalues of this matrix, one has to solve 
the secular equation 

det | A (G t ) — rf | = 
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or 



— X 


1 


1 


1 


1 


— X 
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1 


— X 



det |A(G 2 ) — xl| =det 

= x l -Gx 2 — Sx — 3 = (x+l)2 (a? — 2x — 3) = 0. 



Whence 



*i= +3, 
4 =-l. 



Thus, three roots have the same magnitude (the triple 
degeneration), and all the roots are real because the ma- 
trix A (G 2 ) is Hermitian, i.e. includes only real elements 
and A (G 2 ) = A (G 2 ). 

To find the eigenvalues of the matrix A (G 2 ) we could 
also use the Sachs theorem: 



P (<? 2 ) = a x k + a,x 3 + a 2 x 2 + a 3 x + a 4 = 



a = 1; S x = 0; a x = 0; a 2 = — 6 (according to Corol- 
lary 5 from the Sachs theorem, see p. 45); a 3 — — 8 
(according to Corollary 6 from the Sachs theorem, see 
p. 45). 




s 4 : 



-o > 




a 4 = 3 (-1) 2 2° + 3 (— l)^ 1 = -3. 

Hence: 

P (G 2 ) = x* - 6x 2 — 8x — 3. 

Now let us take account of symmetry. The MOs of a 
tetrahedral molecule are transformed according to irre- 
ducible representations (IR) A x and 2% of the group T d , 
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(For the sake of simplicity we assume that each ligand is 
represented hy only one orbital of s type.) Then in a sys- 
tem consisting of four a-ligands there will be one nonde- 
generated energy level with IR A x and xi = 3, and one 
triply degenerated level with x 2 , 3, 4 — - 1 and IR T~: 

h ( x 2. 3,4= — 1)> 

— ai(^i=+3). 

Thus, knowing the topology and symmetry of a system, 
one can determine the most probable ordering of molecu- 
lar-orbital energy levels and the degeneration structure 
of the energy spectrum. 

Topologically equivalent groups. How is the symmetry 
of a compound related to its topology? Earlier we noted 
that the geometry (symmetry) of molecules may differ 
even when the topology is the same. This statement, how- 
ever, has to be made more precise. To some extent topo- 
logical matrices convey the symmetry of molecules be- 
cause the matrix A, as a rule, commutes with the matrices 
of those ligand permutations that correspond to opera- 
tions of the symmetry point group (R) of the molecule* . 
Consider, for example, a V-shaped molecule of symmetry 
C 2V (say, a water molecule) whose topological matrix 
has the form 

/0 1 0\ 

A = 101 jAs 

\o 1 0/ 

Each symmetry operation of the group C 2v (the symme- 
try group C 2 for a graph) is associated with a certain per- 
mutation of vertices in the topological graph of the mole- 
cule (Table 2). 

The matrices of these permutations 

/0 1\ /l 0\ 

r(P 13 ) = (0 10 and T(/)= 1 
\i 0/ \0 1, 

commute with the matrix A which the reader himself can 
easily demonstrate. 

* The matrices A and B are called commutating if the follow- 
ing relation holds [AB] = AB — BA = 0, i.e. their product is 
independent of the sequence of cofactors. 



* Table 2 



Permutation of graph 
vertices 

Symmetry operation of 



U2 3/ 




E 



group C 2 



The matrices of permutations corresponding to opera- 
tions of molecular symmetry groups also form a group. Let 
us denote it by P n . The groups P n and R are isomorphic. 
In our example the group C 2 is isomorphic to the group P%. 
Since the topological matrix A commutes with the permu- 
tations of the group P n , it will also commute with the 
matrices of transformations of the molecular symmetry 
group. The same group P n will correspond to isomorphic 
symmetry groups. Let us use here the concept of topological- 
ly equivalent groups introduced in 1968 by the German 
chemist H. H. Schmidtke. If two symmetry groups are 
isomorphic and the geometric bodies corresponding to 
them can be obtained from one another through continu- 
ous deformation, such groups are, according to Schmidtke, 
topologically equivalent. If no distinction is made be- 
tween topologically equivalent groups, the topology and 
symmetry of a molecular graph become related in a one- 
to-one manner. Thus, for example, the molecules shown 



in Fig. 35 correspond to the topologically equivalent 
groups Z) 4 and D 2d . The topologically equivalent groups 
possess similar (as to the number, sequence, and degen- 
eration structure) spectra of molecular-orbital energies 
(Table 3). 

Structural matrix. Strictly speaking, the foregoing argu- 
ments are correct only if each atom provides one orbital to 
form a chemical bond. When this condition is not satisfied 




Fig. 35 
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the matrix A, as was shown by Schmidtke, should be 
replaced by a matrix of higher dimension which he called 
a structural matrix. Then the group P n will represent 



Table 3 



Eigenvalues of topological 
matrix A 


Irreducible representations of groups 


Dk 
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Ai 


X 2i 3 = 

x 4 =— 2 
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E 
B t 



the permutations of AOs. For example, the structural 
matrix N corresponding to the tetrahedral system L 4 




Fig. 36 



whose each vertex accommodates three p-AO (Fig. 36) 
has the form 



N= 
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Here the dimensionless parameter x allows, for only 
the n interactions between ligands. The spectrum of eigen- 
values of the matrix N and the corresponding IR groups 
T d are given in Table 4. 



Table 4 



Eigenvalues of matrix N 


IR groups 
T d 


*l = 4 + X 
x 3 = — 1 — x 

*4. 8 = - 4- ±Y t" - 3x + x2 


Ai 

E 
Ti 



From the data of Table 4 we obtain the following se- 
quence of MO energies: 

e(a x ) < e (t 2 ) < e (e)< e (fj < e (f 2 ). 

The calculations of tetrahedral systems (B 4 H 4 , B^J 2 
etc.) confirm the result obtained. 

Note that although the matrix N is of rather high 
dimension, it is more convenient to use it to get the 
qualitative and semiquantitative information than to 
calculate the energy levels directly by the quantum- 
mechanical methods because the matrix N is substan- 
tially "rarefied" with zero elements. 

Topology of Molecular Charge Distributions 

The development of chemistry has both led to, and been 
made possible by, the evolution of certain primary con- 
cepts. These concepts, without which there would be 
neither correlation nor prediction of the observations of 
descriptive chemistry, are: (1) the existence of atoms of 
functional groupings of atoms in molecules as evidenced 
by characteristic sets of properties; (2) the concept of 
bonding; and (3) the associated concepts of molecular 
structure and molecular shape. These concepts logically 
(but not historically) are consequences of fundamental 
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topological properties of the charge distribution (electron- 
ic and nuclear) in a molecular system. In terms of the 
Born-Oppenheimer approximation the electronic distri- 

bution p (r) is the scalar field defined in the real three- 
dimensional space with Euclidean metric. 

The universal topological properties of p (r) are charac- 

terized by its gradient field Vp (r). (See Fig. 37 showing 

the field p (r) and the trajectories of its gradient in the 
B 2 H 4 plane of diborane.) 

Critical points. The properties of this field, and hence 
the principal characteristics of a charge distribution, are 
totally determined by the number and nature of its critical 
points, i.e. the points at which the field vanishes. Exactly, 

in the neighbourhood of the critical point r c 
Vp(r)U -* =0. 

I r=r c 

K. Collard and G. Hall have demonstrated the utility 
of orthogonal trajectories, i.e. the paths traced by the 
gradient vectors of a scalar field. The gradient path direc- 

tion is taken to coincide with that of growing p (r) values. 
Only a single gradient path traverses any given point in 
space. All gradient paths are orthogonal to the isodensity 
(contour) lines. Also, all the paths originate and termi- 
nate at the critical points. Besides, Collard and Hall 
indicated that the properties of a critical point can be 
uniquely defined by two numbers, labelled rank and signa- 
ture. R. F. W. Bader and coworkers have adopted them 
in their classification of charge distributions. 
A bit of differential geometry... The set of nine second 

derivatives of p (r) 

forms a so-called Hessian matrix A (named after the Ger- 
man mathematician Ludwig Otto Hesse (1811-1874)) of 

p (r) at the critical point. The properties of this point are 
determined by diagonalizing the Hessian matrix in the 

neighbourhood of r c . The rank of the critical point is equal 
to the number of nonzero eigenvalues of A. Usually the 
charge distributions in molecules possessing the equilibri- 
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Fig. 37 



um nuclear configurations are characterized by the criti- 
cal points of rank 3. The three eigenvalues determine the 
three orthogonal gradient paths that either originate or 
terminate at the critical point. 
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The signature is the number of excess positive eigen- 
values over negative ones. 

Types of critical points. Graphically, the critical point 
is characterized by its phase portrait, i.e. the pattern 
of trajectories traced out by the gradient vectors in its 




Fig. 38 



neighbourhood. For critical points of rank three, four 
characteristic phase portraits are found, corresponding to 
the four possible signatures, — 3, — 1, +1, +3. Each of 
the four possibilities labelled as [rank, signature], e.g 
[3, — 1], denotes the presence of a particular element of 
molecular sti ucture. 

A [3, — 3J critical point. In this case all three eigen- 
values are negative and hence there is a local maximum in 

— *■ 

p (r). All the gradient paths in the neighbourhood of a 
[3, —3] critical point terminate at that point (see Fig. 38a). 



8-01102 
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The maximum in p (r) at the position of a nucleus is not 

a true critical point as Vp (r) does not vanish. However, 

near the nuclear cusp p (r) behaves topologically as if it 
were a true [3, — 3] critical point (see also Fig. 37). 

A (3, —1] critical point. In this case p (r) has a saddle 
point located between neighbouring atoms. All trajec- 
tories forming a surface terminate at the critical point 
(Fig. 386). On this surface, the [3, — 1] point has the 

properties of a two-dimensional maximum in p (r), i.e. it 
behaves as a [2, —2] critical point. The eigenvector cor- 
responding to the positive eigenvalue defines a unique 
axis for the [3, —1] point, which is perpendicular to the 
surface at the saddle point. The two gradient paths defin- 
ing this axis originate at the critical point. These paths 

also define a line along which p (r) is maximum with respect 
to the directions perpendicular to that line. Such a line 
is called a bond path. 

It has been demonstrated that the network of bond 
paths defined by the unique axis of the [3, — 1] critical 
points coincides with the network formed by linking to- 
gether those pairs of atoms which are presumed to be 
bonded chemically. 

Two remaining critical points... These points, [3, +1] 
and [3, +3], appear in actual molecular graphs. (Accord- 
ing to R. Bader a molecular graph is a network of bond 
paths found for a given molecule.) If three or more nuclei 
form a ring of bonded atoms, then a [3, +1] point is to 

be found within the ring's perimeter. The associated ring 

— +■ 

surface is one on which p (r) possesses a minimum value 
at the critical point [3, +1]; the gradient paths forming 
that surface originate at the saddle point [3, +1] (Fig. 38c). 
Figure 38d illustrates the phase portrait of an atom with 
a single neighbour (all trajectories terminate at the 
[3, — 3] point). 

If four or more nuclei form a bonded cage, a true three- 
dimensional minimum in p (r) must exist inside that 
cage; all gradient paths within the cage originate at such 
a critical point. 

The number and type of critical points of rank three 
which can coexist in a molecular system is governed by the 
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following expression (the partial case of the Poincare- 
Hopf relationship): 

n — N + r — c = 1, 

where re is the number of nuclei, i.e. the number of pseu- 
do-[3, —3] critical points, N is the number of bonds, i.e. 
the number of [3, — 1] critical points, r is the number of 
rings, i.e. of the [3, +1] critical points, and c is the 
number of bonded cages, i.e. of the [3, +3] critical points. 
For example, the (re, N, r, c) index characterizing the 
structure of P 4 is (4, 6, 4, 1) and of diborane (8, 8, 1, 0). 

Clearly, the set of numbers (re, N, r, c) is a useful index 
of the major topological properties of a molecular charge 
distribution. 

When structure is changed... The making and breaking 
of a chemical bond is topologically a discontinuous process, 
and the associated change in structure is therefore abrupt. 
It has been pointed out that the analysis of the discontin- 
uous change in the topological characteristics of a molec- 
ular charge distribution resulting from the discontinuous 
change in the nuclear coordinates is given by the catastro- 
phe theory of Rene Thorn. However, the discussion of 
that theory lies outside the scope of our book. 



Chapter Seven Live Pages 

of the Past 



In this chapter we shall tell of the origination of the 
graph theory and its first utilization for solving various 
chemical problems. At the same time we shall consider 
some present-day problems pertaining to chemical topol- 
ogy and having, as will be shown below, a long history. 

Analysis Situs and Origination 
of Graph Theory 

Walk over Konigsberg bridges. The need for the devel- 
opment of new «nonmetric» geometry was voiced first 
back in the 17th century. Thus, the great German mathe- 
matician and philosopher G. W. Leibniz wrote on 
September 8, 1679 in a letter to the Dutch physicist C. Huy- 
gens: "I am chagrinned with algebra because it provides 
neither short nor best ways to construct geometry, and 
I think therefore that we need still another analysis, 
properly geometrical ... permitting a position to be expres- 
sed just as directly as algebra expresses a quantity". 

The historians disagree on what was actually meant 
by Leibniz, whether it was what later became topology, 
or the branch of mathematics called vector analysis. 
Everybody is unanimous, however, concerning one thing: 
the origination of topology and the graph theory is asso- 
ciated with the work of the great mathematician of the 
18th century Leonard Euler (1707-1783), a Swiss who lived 
in Russia for a long time and who made significant contri- 
butions to various branches of pure and applied mathe- 
matics. During his lifetime, half of which he was absolute- 
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ly blind, „Euler published 530 books and articles and 
after his death, the remaining manuscripts were being 
published for almost 50 years more, so that the total num- 
ber of his works approached 900. 

The article that is of interest to us now was submitted 
by Euler to the St. Petersburg Academy of Sciences on 
August 26, 1735. Wrote Euler: "In addition to that part of 
geometry where magnitudes are treated and which was 
zealously elaborated at all times, there is another, quite 
unknown so far, part of geometry to which Leibniz allud- 
ed for the first time, calling it geometry of position.... 
I believe it will be useful to demonstrate here as an 
example of geometry of position a method I found to 

solve problems of that kind One of such problems 

which is well known boils down to the following. 

In the city of Konigsberg, East Prussia, there is an 
island named Kneiphoff . It is surrounded by a bifurcating 
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river spanned by seven bridges (Fig. 39). The question 
is whether one can walk over all bridges crossing each of 
them only once." 

Euler proved such a stroll over Konigsberg bridges 
impossible. Moreover, he generalized his conclusion, hav- 
ing demonstrated that the required itinerary exists only 
for the connected graph each of whose vertices is incident 
to an even number of edges. As to the graph representing 
the problem of Konigsberg bridges (Fig. 40), it is con- 
nected, but its vertices are of only odd degrees. 

Euler proved another topological theorem (1750, pub- 
lished in 1758) named after him: for every convex poly- 
hedron homeomorphic to a sphere, 

$ — & + & = 2, 
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where $i is the number of vertices, & the number of 
edges, and 3" the number of faces. 

To some of our readers this formula may remind Gibbs' 
phase rule 

/ -ft + r = 2, 

where / is the variance of the system, i.e. the number of 
independent variables specifying its state (pressure, tern- 




Fig. 40 

perature, component concentrations of all phases); k is 
the number of components, and r is the number of phases 
of the system. 

The phase rule plays an important part in chemistry. 
For example, it helps to answer the controversial question 
of long standing: whether the solubility of solids in 
liquids depends only on temperature, or also on pressure. 
Inasmuch as such a system is composed of two phases 
(liquid solution and a solid) and two components (solvent 
and solute), then 

f = 2 + 2 — 2 = 2. 

Thus, the considered system is divariant, and the sol- 
ubility depends on two variables, temperature and 
pressure. 

Graphs and electricity. After Euler the graph theory 
was rediscovered several times. In 1847 G. Kirchhoff devel- 
oped the theory of trees to solve the problem about cur- 
rents flowing in each conductor and in each loop of an 
electric circuit. Kirchhoff replaced the electric circuit 
by the graph and developed the procedure to solve the 
system of equations defining the current intensity, analys- 
ing only the so-called skeleton trees of the initial graph. 
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A few words on terminology. In the 19th century the 
great contributors to the development of the graph theory 
were many eminent mathematicians: A. Cayley (see 
below^; K. Jordan, J. J. Sylvester and others. At the 
same time, J. B. Listing, F. Klein, A. Mobius, B. Rie- 
mann, A. Poincare and others developed topological 
methods in their works. The term "topology" itself was 
introduced in 1836 by J. Listing, but was put into 
practical use only in the late 20s of the 20th century; 
up to that time the Latin term analysis situs was prefer- 
red ("analysis of position"), which is traced back to 
Leibniz. 

However, before mathematical methods could be em- 
ployed in chemistry, the latter had to reach sufficient theo- 
retical maturity. Thus, the graph theory and some other 
divisions of discrete mathematics could be applied to 
the solution of chemical problems only after the chemical 
structure; theory had emerged. In the next section we 
shall briefly discuss some important landmarks in the 
development of ideas concerning the structure of chemical 
compounds. 

From Ideograms to Structural Formulas 

Symbol of substance. Since ancient times people have 
used various symbols for conventional description of real 
objects and phenomena. On the walls and columns of old 
Egyptian temples one can see the signs of water, gold, 
silver, etc. Figure 41 shows the signs that were used by 
three ancient civilizations, Chinese, Egyptian, and Greek, 
to denote water and alum. These signs are of different 
nature; sometimes they "picture" a substance or some of its 
properties (see, e.g., the Egyptian sign for water), whereas 
in other cases ideograms convey some secret meaning 
which is often very difficult to decipher. 

The ideographic signs were used by medieval alchemists 
although sometimes they resorted to letter abbreviations 
(Latin, as a rule) of substances. The medieval nomencla- 
ture was very jumbled. At the end of the 18th century the 
Swedish chemist T. 0. Bergman (1735-1784) and the 
French scientists A. F. Fourcroy, L. B. Guyton de Morveu, 
and A. L. Lavoisier proposed to put it in order. However, 
for some reasons, that we are not able to discuss now, 



their system had not become popular. The descriptive 
signs of some compounds and the chemical reaction sym- 
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bols used by the French chemists were very complex and 
ambiguous. 

An interesting version of graphical notation for chemi- 
cal formulas was proposed by the English scientist 
W. Higgins. Figure 42 illustrates the formulas for some 
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nitrogen oxides, taken from his treatise A Comparative 
View of the Phlogistic and Anti-Phlogistic Theories pub- 
lished in London in 1789. The symbol P denotes nitrogen, 
and the symbols a, b, c, etc., oxygen atoms. The numbers 
indicate the atomic interaction intensities. To each atom 
Higgins assigned three conditional "force units". Then the 
total force of atomic interaction equals 3n, where n is 
the number of atoms in the molecule, and the value 
Znlk is alloted to each bond (k is the number of bonds). 
It should be noted that Higgins had not defined the term 
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"bond", it> was introduced into science only three quarters 
of a century later. 

Emergence of chemical atomism. The beginning of the 
19th century was marked by the emergence of atomic 
theory developed by the famous English scientist John 
Dalton (1766-1844). According to Dalton, atoms combine 
in multiple proportions. Thus, if two atoms form only one 
compound, they associate in the 1 : 1 ratio; if two atoms 
form two compounds, their ratios in these compounds 
are respectively 1 : 1 and 1 : 2, and so on. Accordingly, 
Dalton adopted the formula HO for water, CO for carbon 
monoxide, C0 2 for carbon dioxide, and so on. From these 
formulas he determined the atomic masses of elements. It 
is clear that frequently both atomic masses and chemical 
formulas of Dalton were not correct. 

Besides, Dalton proposed new symbols for chemical 
elements, as well as a new way of graphic representation 
of compounds. Some of the symbols and formulas taken 
from his treatise A New System of Chemical Philosophy 
(1808) are shown in Fig. 43. Note that Dalton sought 



Q - oxygen 



aluminium 



^ — sulphur 



— potassium 



^) — hydrogen 




— alum 



Fig. 43 



to present his graphical formulas in the most ordinary 
and symmetric fashion. 
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Later the Dalton method of determining the 'formulas 
of chemical compounds was castigated by the greatest 
Swedish chemist J. J. Berzelius (1779-1848)*. 

However, the method of determining the element compo- 
sition of substance, proposed by Berzelius, also failed 
because the reliable methods of determining molecular 
masses had not been elaborated yet. The existing methods, 
as a rule, led to confusion resulting in its turn in subse- 
quent numerous errors in determining the molecular 
masses of the whole classes of organic compounds. 

The first adequate method of determining molecular 
masses was the method based on the measurements of 
vapour density of a substance. That method was a direct 
consequence of the Avogadro molecular theory developed 
in 1811-1821. However, the ideas of the Italian scientist 
were recognized only after the congress of chemists in 
Karlsruhe (1860), primarily due to the efforts of his com- 
patriot S. Cannizzaro. The latter completed reforming 
Dalton's atomic theory and opened the way to the 
emergence of the chemical structure theory. That theory, 
as well as many significant concepts of structural chemis- 
try (e.g. valence, or chemical bonding), could appear only 
when an end was put to the confusion and arbitrariness 
in derivations of chemical formulas and determinations 
of molecular masses. But the structural theories which 
appeared before the 60s of the last century were still 
marked with the imperfection of the earlier atomic 
concepts. 

Variation on four themes. The most significant of these 
theories is the theory of types proposed by the eminent 
French chemist Ch. Gerhardt (1816-1856) in 1853. Accord- 
ing to Gerhardt, each organic compound belongs to one 
of four "types": H 2 , HC1, H 2 0, and NH 3 . Replacing 
hydrogen in them by various radicals, one can obtain 
hydrocarbons (H 2 type), alcohols, ethers, acids (H 2 
type), etc. On the basis of his theory Gerhardt developed 
a harmonious systematization of organic compounds^using 
extensively the concept of homology; he also predicted, 
and at a later date synthesized, a number of anhydrides 
of monobasic acids. However, as the German chemist 



* Besides, in 1813-14 Berzelius proposed the letter designations 
of elements which on the whole are used even now. 
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H. Kolbe aptly noted, "nature does not limit itself with 
performing variations on four themes: water, hydrogen, 
hydrogen chloride, and ammonia". In fact, the number of 
types had to be later considerably increased, so that the 
theory had lost much of its initial simplicity and clear- 
ness. Moreover, it could not cover many facts related to 
the addition and elimination reactions and to isomerism. 

Kekule's formulas. In 1857 the German chemist A. Ke- 
kule deduced that oxygen was bivalent, nitrogen and 
phosphorus trivalent, and, what is most important, car- 
bon tetravalent. He regarded valence as the number of 




Fig. 44 



"affinity units" capable of mutual saturation when atoms 
combine into a molecule. Similar ideas were developed 
by H. Kolbe (1875) and the English chemist E. Frank- 
land (1852). 

Figure 44 is a graphic representation of methyl ether 
of formic acid according to Kekule (1859). 

Referring to his formulas, Kekule warned that they do 
not express either the spatial distribution of atoms or 
any metric properties of molecules; the formulas only 
described the "relationship of binding atoms". 

Three aspects of chemical structure theory. By the 
middle of the 19th century the time was ripe for radical 
changes in organic chemistry since none of the contempo- 
rary theories could generalize the tremendous amount of 
accumulated experimental data. And the new theory was 
not long in appearing. It appeared in 1861; that was the 
theory of chemical structure produced by A. M. Butlerov. 

The term "chemical structure" was defined by the Russi- 
an scientist in different ways; first, as "chemical" bonding 
or mode of mutual binding of atoms in a complex body", 
later, as "distribution of action of affinity a certain 
amount of which is possessed by each chemical atom", and 
still later, as "mode of mutual chemical bonding of ele- 
mentary atoms in a molecule". In their essence, all these 
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definitions are quite similar. In modern scientific terms, 
Butlerov represented every chemical compound as a topo- 
logical graph, having associated the properties of a mole- 
cule, its "chemical nature", with "the nature of elementary 
constituents, their number and chemical structure". 
Butlerov especially emphasized the thesis about recipro- 
cal influence of atoms in a molecule. 

An important aspect in the development of the chemical 
structure theory was its extension to unsaturated com- 
pounds. In 1870 Butlerov demonstrated that the concept of 
multiple bonding formulated earlier by A. Crum Brown 
and E. Erlenmeyer and furnished an adequate interpreta- 
tion to all cases of isomerism, known at that time, in the 
unsaturated compound series. The introduction of the 
multiple bonding concept into structural chemistry made 
it possible in many cases to solve the problem of distri- 
bution of "the power of chemical affinity" among atoms. 

Thus, we can distinguish at least three aspects of the 
chemical structure theory: 

(1) topological aspect, when neither the multiplicity 
of chemical bonding, nor any of its energy or geometric 
characteristics, and reciprocal influence of atoms are 
considered; only the topological graph of a molecule is 
specified; 

(2) mixed valence-and-topological aspect, when the no- 
tion of chemical bonding is generally characterized by 
indicating its multiplicity, taking into account formal 
valencies of atoms; molecules are thus represented as 
graphs with multiple edges; 

(3) proper chemical aspect, when individual character- 
istics of atoms in a molecule, their reciprocal influence, 
specific nature of bonds, etc. are taken into consideration. 

Riddle of isomerism. The chemical structure theory 
made it possible to solve many urgent problems of chem- 
istry of that time, in particular the problem of isomerism. 
The vestiges of ideas about isomerism can be found in 
the atomist theories of philosophers of Ancient Greece. 
Many years later, in the 18th and the beginning of the 
19th century, some scientists (J. Jungius, J. Dalton, 
A. von Humboldt) advanced the theory of the existence 
of substances of identical composition, but possessing 
different properties due to different arrangement of atoms 
in their molecules. In 1811 the French physicist and chem- 
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ist J. Gay-Lussac proved experimentally the existence 
of isomers. However, the isomerism concept was introduced 
into chemistry only in 1830s by Berzelius who defined 
isomers (from Greek ioo£, equal, and fiepo£, a part, share) 
as "substances possessing identical chemical composition 
and molecular weight, but different properties". 

The underlying cause of isomerism was disclosed only 
after the structure theory was produced. Owing to the 
work of Butlerov and his followers the isomers of various 
classes of organic compounds, alkanes, olefins, alcohols, 
etc., were investigated. In 1867 the German chemist 
W. Korner, a student of A. Kekule, introduced the terms 
ortho-, meta-, and para-isomer for disubstituted deriva- 
tives of benzene. 

By the middle of 1870s the problem of the number of 
isomers for a compound of a given composition acquired 
more significance, having engaged the attention of both 
chemists and mathematicians. 

Mathematics of Structural Formulas 

Pioneer. One of the first mathematicians who paid 
attention to the chemical structure theory was the 
outstanding English scientist Arthur Gayley. In 1875 he 
employed the graph-theoretical methods for solving the 
problem of counting isomers of C n H 2n+2 alkanes with the 
given n. Initially, Gayley tried to derive the general 
formula for determining the number of such isomers each 
of which he represented as a topological graph (a tree, to 
be more exact). These attempts, however, failed. Then 
using the results of his earlier 1857-59 investigations, 
Gayley developed the method of counting so-called rooted 
trees with n vertices. A tree is called rooted if one of its 
vertices is marked, for example, somehow coloured. 
Figure 45 shows nine rooted trees with five vertices 
(n = 5)* . Then Cayley divided the set of rooted trees into 
two groups, centric and bicentric trees. 

The centre of a graph is such a vertex v whose maximum 
distance from any other vertex is the least of all possible 
distances. Figure 46 illustrates centric and bicentric trees. 

* To satisfy the requirement for carbon to be tetravalent, Cayley 
considered only the trees whose vertices have the degrees not ex- 
ceeding four. 
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Let us get back to our example (n = 5). Ampng nine 
trees pictured in Fig. 45 there are five centric trees, a, d, 
e, g, h (the centre of the graph may not coincide with the 
coloured vertex), and four hicentric ones, b, c, /, i. 




'a) (bi (c) (d) (e) 




(f) (g) (h) (i). 



Fig. 45 

It is easy to see that all bicentric trees are isomorphic 
(the fact that one of the vertices is coloured is not taken 




ceniric bloenirla 
Fig. 46 



into account). As to the five centric trees, a, d, h and 
e, g are isomorphic within each group. 

Thus, the set of nine rooted trees with five vertices can 
be subdivided into two subsets of nonisomorphic centric 
trees and one subset of bicentric trees. According to Cay- 
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ley, the .total number of isomers of saturated C H 
hydrocarbons is equal to the total number of nonisomor- 
phic rooted trees, i.e. 3 (=2 + 1). 



Table 5 
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15 
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10 
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12 
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357 
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13 


419 


380 


799 


802 



In Table 5 Cayley's results for the first 13 members of 
the homological series of alkanes are compared to the 
present-day data. As is seen from the table, Cayley's 
method provides erroneous results beginning with n — 12. 

Beachhead. However, basically this discrepance does not 
disparage the explorations of the English mathematician. 
Cayley was the first to represent a molecule as a topologi- 
cal graph and to start developing the methods of isomer 
enumeration on the basis of the graph theory. Moreover, 
Cayley forwarded the notion of enumerative polynomial 
for rooted trees, i.e. the polynomial whose coefficients A t 
specify the number of rooted trees with i vertices: 

1 -f A x x + A 2 x 2 + . . . 

The abstract of Cayley's paper devoted to enumerating 
alkane isomers was submitted to the German Chemical 
society (Deutschen chemischen Gesellschaft) by the emi- 
nent German chemist C. Schorlemmer and was published 
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in Berichte der Deutschen chemischen Gesellsclfaft; the 
paper attracted the attention of both chemists and mathe- 
maticians. 

Later, beginning with 1930s analytical methods of 
finding the coefficients A t were elaborated*. But, as it 
was noted by the greatest expert in the graph theory 
F. Harary, "the methods of enumeration in combinatorial 
analysis may be regarded rather as art than as science". 

A. Cayley was also the creator of the matrix theory; 
he made an essential contribution to the development of 
the group theory, i.e. those branches of mathematics which 
were later extensively used in physics and mathematics. 
Moreover, Cayley was the first to indicate the relation- 
ship between the point groups of symmetry and the 
permutation groups (see Chapter 6). 

Cayley's theorem states: 

Any finite group G of order n is isomorphic to a certain 
subgroup of the permutation group S n . 

Note that the considered group G of the re-th order is 
isomorphic not to the group S n itself but to its subgroup, 
since S n contains more elements than G does. For example, 
the symmetry group of an equilateral triangle C 3v includes 
six symmetry operations (E, C 3 , C\, 3cr ) whereas the 
group iS s includes 6! permutations. Consequently, C 3D 
and S 6 are not isomorphic, but the groups C 3o and S 3 are. 

Since 1928 permutation groups have been extensively 
used in quantum mechanics of poly-electronic systems 
(atoms, molecules, solids) and during the last decade 
they, together with the graph theory, found application 
in the theoretical treatment of some intramolecular rear- 
rangements. Due to the significance and obvious usefulness 
of the results obtained thereby, we shall discontinue for 
a while the chronological sequence of presentation and 
tell of certain methods of description of spatially nonrigid 
compounds. 

Ligands in Motion 

The study of some nonrigid molecules carried out 
during the last 20 years or so revealed various intramole- 

* The most significant papers were published by J. H. Redfield 
(1927) and G. Polya (1937); also see the monograph Graphical 
Enumeration by F. Harary and E. M. Palmer, Academic Press, 
N. Y., London, 1973. 
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cular rearrangements of ligands taking place in such 
molecules and resulting in the ligand interchange which 
occurs without a break of their bonding with the central 



£ 




Reaction coordinate 



Fig. 47 

atom. Various mechanisms of such a ligand interchange 
were suggested. An example is provided by so-called... 

Berry's pseudorotation. In the case of a trigonal bipyr- 
ramid molecule, e.g. PF 5 , this mechanism consists in an 




(b) 

Fig. 48 




interchange of two equatorial ligands with two apical 
ones, with the third equatorial ligand keeping its position. 
Due to such a rearrangement the molecule rotates in 



9-01102 
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pace, hence the term "pseudorotation". Figure 47 shows 
energy variation vs the reaction coordinate in the process 
of Berry's pseudorotation. The peak in the curve corre- 
sponds to the transition state with tetragonal pyramid 
configuration. In many penta-coordination compounds 
the energy barrier is low and can be overcome at higher 
temperatures. 

The pseudorotation phenomenon is vital to structural 
chemistry of phosphorus, silicon, transitional metals, 
and some other elements. Apart from pseudorotation, 
there are other rearrangement mechanisms, for example, 
trigonal twist (Fig. 48a) examined in dynamic stereochem- 
istry of icosahedral carboranes, "tourniquet" mechanism 
(Fig. 48ft), etc. 

The theoretical investigation of these reactions re- 
quires quantum-mechanical methods, in particular, the 
study of chemical bonds in initial, final, and intermediate 
compounds, as well as the consideration of nuclear mo- 
tions. Yet frequently the important information can be 
obtained without analysis of electronic structure of 
molecules and investigation of actual motion of nuclei, 
only resorting to the graph theory and employing the 
group-theoretical conceptions. Here we should define two 
terms: permutation isomers and permutation isomerism 
reaction. 

Permutation isomers are identified as the com- 
pounds having the same qualitative and quantitative 
composition, the identical set of ligands, the identical 
molecular skeleton, and differing only by the dis- 
tributions of ligands over skeletal positions in the 
molecule* . 

Let us label the ligands: {l ly l 2 , • • ., l n } and the skele- 
tal positions in the molecule: {s^ s 2 > • • •> s n}- Then 
each permutation isomer will have a corresponding (2 X 
X n) matrix: 

fl 2 3 4 5 ... n\ 

\s ) \i j k t u ... m)' 

The first row of the matrix consists of ligand indices 
(numbers) written in the growing order, and the second 

* The term "skeletal position" refers to the "place" in the vicinity 
of the central atom at which a ligand can be found. 
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one of indices (numbers) of skeletal positions. As an 
example, let us consider a trigonal bipyramid system 
with a definite enumeration of skeletal positions: 




Denote five ligands as follows: i x as A, Z 2 as B, l 3 as C, 
1 4 as D, and Z 5 as E. Then the permutation isomer 




will be represented by the matrix 



[s )~\5 2 1 3 AJ 



implying that the first ligand (.4) is in the fifth position, 
the second (B) in the second, the third (C) in the first and 
so on. The isomer of the type 

E 

is adopted as the standard one. 

When one permutation isomer of a certain molecule 
turns into another, differing from the first one, we say that 
the permutation isomerism reaction takes place. 

The following rearrangement may serve as an example of 
such a reaction: 




9* 
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The symbol p t = (123) (4) (5) denotes a rearrangement 
that brings the ligand that resided initially at position 1 to 
position 2, the ligand that resided initially at position 2 
to position 3, and the one that was at position 3 to posi- 
tion 1. Such a cyclic permutation is designated by the 
symbol (123). This expression is referred to as a cycle. 
The length of a cycle, i.e. the number of elements in it, 
is equal to 3 here. Cycles of unitary length imply that 
the ligands residing at positions indicated do not change 
these positions as the reaction proceeds. In our case the 
ligands at positions 4 and 5 do not shift anywhere (cycles 
(4) and (5)). It should be noted that the arrows pointing 
to the ligand displacement direction are conditional and 
may not correspond to the actual rearrangement mecha- 
nism, the latter can be determined either experimentally, 
or theoretically, by quantum-mechanical methods. 

Why weaving a tangled web? the reader may wonder. It 
is well-known that n objects can be permutated in w! 
ways. Is it not simpler to write out all possible permu- 
tations of ligands (no trouble if n is large: a computer 
will help out) and then to distinguish those with the 
lowest activation barriers? 

This can be done only if a molecule has no symmetry 
elements. But in most cases the number of different reac- 
tions of permutation isomerism is less than ra! . And there 
are definite reasons for that. 

First, some permutations correspond to the rotation of 
a molecule about one of its symmetry axes, for example, 

-:■ : * A 

r,.=(/)(234)(5)! 
B ^> D 




and do not lead to permutation isomerism reactions be- 
cause in these cases a molecule just turns in space. 

The set of all permutations involving the series of 
indices of skeletal positions and corresponding to the 
rotation of a molecule as a whole forms a group (denoted 
usually by R) which is actually a subgroup of group S n 
of all permutations of n ligands. 
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Any symmetry operation (rotation as well as reflection) 
can be associated with a certain permutation. In other 
words, the symmetry operations of the molecular point 
group G generate a permutation group of skeletal position 
indices (G s ). Note that G s is also as a rule nonisomor- 
phic to the group S n (see above). 

Second, if all ligands are identical in their nature and 
differ only in their numbers, all possible isomerism reac- 
tions prove to be indistinguishable due to the symmetry 
of a molecule. 

Nevertheless such permutations can be distinguished, 
for example, using the method of nuclear magnetic reso- 
nance (NMR), e.g. the geometry of a phosphorus penta- 
fluoride molecule was found by the electron diffraction 
method and from the analysis of IR spectra: 



However, the detailed studies of NMR spectra led to 
the conclusion about the interchange of apical and equa- 
torial fluorine atoms. Since the apical and equatorial bonds 
P— F in PF 5 have different lengths, the NMR spectrum 
can be expected to consist of two doublets with integral 
intensity ratio 2 : 3. Actually, the observed NMR spec- 
trum of PF 5 (at 19 F) has one doublet caused by spin-spin 
interaction of P and F nuclei, i.e. the splitting of signals 
from axial and equatorial fluorine atoms is absent. 

Thus, far from all permutations can be treated as permu- 
tation isomerism reactions. The distinguishable permuta- 
tions should be selected. The criterion of indistinguish- 
able permutations of ligands was proposed in 1971 by the 
American chemist W. G. Klemperer; it can be stated 
as follows: 

Two permutation isomerism reactions p n and p m 
are indistinguishable (identical) in the case of 
completely symmetrical (achiral) surrounding of 
the central atom if for a certain permutation g t 6 G s 
the equality 




p r ( symmetry D 3h ) 



F 



Pn = gl l Pmgi 



(36) 



holds, 



Otherwise the reactions p n and p m can„be distin- 
guished. 

What is chirality? The term chirality denotes the 
property of a molecule to be nonidentical to its mirror 
image. All coordination polyhedrons with identical li- 
gands are achiral. When the ligands are different, the chiral 
compounds may be formed. Thus, a tetrahedron with four 
different ligands is chiral. A trigonal bipyramid and 
tetragonal pyramid may be (but not necessarily will be) 
chiral if three (or more) ligands are of different types, 
e.g. 



L L 




Chiral achiral 



Subsequently we shall confine ourselves to the case of 
intramolecular isomerism in achiral molecules. 

Figure 49 illustrates the relationship between two in- 
distinguishable (according to Klemperer) permutation 
isomerism reactions [g t = (1) (234) (5) represents the 
rotation about the axis C 3 passing through the ligands AE\ ■ 

Counting indistinguishable reactions. The number of 
indistinguishable reactions is determined by means of the 
enumeration polynomial whose coefficients are calculated 
by the graph-theoretical methods and the methods of per- 
mutation group theory. We cannot describe here the pro- 
cedure of finding the coefficients of the enumeration poly- 
nomial, and so we shall consider only several concepts that 
are necessary for translation of mathematical results into 
the language of chemistry. 

The set of permutations {p t } satisfying the condition 
(36) forms the coset of the group S n (and its subgroup G s ). 
The total number of such cosets minus the number of sets 
consisting of the elements of group R is equal to the num- 
ber of formally distinguishable permutation isomerism 
reactions for an achiral molecule. 

It can be demonstrated that permutations relating to 
the same class have the same cycle structure, or expressed 
^otherwise, are of the same cycle type. The following 



example 'may serve as an illustration: two symmetry 
operations'C 3 r of group D 3h correspond to two permutations 
of the same class: 



(1) (2 3 4) (5) 
(1) (2 4 3) (5) 



Both of them contain two cycles of unitary length and 
one cycle of length 3 each. The cycle structure of these 
permutations is expressed by the symbol (2, 0, 1, 0, 0), 
where J the ordinal number of the digit in parentheses 
coincides with the cycle length and its magnitude with the 
number of cycles of a given length. In the symbol (2, 0, 
1, 0, 0) the first place is occupied by digit two because 
both permutations contain two cycles of unitary length 
each and so on. In the general case the cycle structure 
symbol of permutations of a given class has the form: 
(/n In • • •> /n)» where / x is the number of cycles of length 1, 
/ 2 is the number of cycles of length 2, etc. Each distin- 
guishable permutation isomerism reaction belongs to its 
specific class and cycle type. 

Now we can write the enumeration polynomial F in the 
general form defining the permutation isomerism reac- 
tions: 

F= 2 j n a»a% . . . ain, (37) 

where A) lt j lt . . j n is the] numb er of jdisting uishable 
reactions of cyclic type / 2 , . . ., /„); thejprime signifies 



A A D n 



p£^?~ B WW)-*', c , 
E 

I 



C 




(MX2)(j)(5)=p„ 



Fig. 49 



that the permutations corresponding to the rotation of the 
molecule as a whole are not taken into account; a t are 
certain conditional variables; summation is carried out 
over all cyclic types. 
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Fig. 50 



Using the graph enumeration methods (so-called P6- 
lya's theorem), Klemperer derived F for the case of 
achiral trigonal bipyramid molecule (D 3h ) in the form 

F = 3a\a\ + 2a\a\ + 2a\a\ + Sa\a\ + 3a\a\ + 2a\ (38) 
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Let us "decipher" this expression. As an example, let 
us consider the first term 3a\a\. The subscript (i) of the 
variable a? denotes the cycle length, the superscript (k) 
the number of cycles of length i. Then the first term in 
the expression (38) means that there are three different 
permutations of cyclic type (3, 1, 0, 0, 0). 

Figure 50 shows all 15 distinguishable permutation 
isomerism reactions for achiral pentacoordination mole- 
cule of symmetry D 3h . 

Klemperer examined the other cases as well: (1) when 
the molecule is chiral, and (2) when the molecule (either 
chiral or achiral) is not fixed in a crystal lattice but per- 
forms rotations and displacements (in fluids). 

Thus, valuable information concerning the possible 
nature of permutation isomerism reactions was obtained 
only by the methods of the graph theory and group theory, 
without analysis of electronic and nuclear densities. 

Cayley's Followers 

The works of A. Cayley on isomer enumeration had 
drawn attention of many chemrsts and mathematicians. 
In 1882 Cayley delivered in< the USA a series of lectures 
on application of mathematics in structural chemistry. 
A few years earlier, in 1878, the English mathematician 
J. J. Sylvester who founded the first mathematical jour- 
nal inj America (Amer. J. Math.) published in it an arti- 
cle dealing with the application of algebraic methods in 
the valence theory. True, Sylvester was not occupied with 
the problem of isomer enumeration. His investigations 
taken up later in 1901-1904, by P. Gordan in Germany 
and V. G. Alekseev in Russia laid a foundation for the 
spin valence theory created a quarter of a century later. 
The authors of that theory (W. Heitler, F. London, 
M. Born, H. Weyl, Yu. B. Rumer, and others) were fa- 
miliar with the investigations of Sylvester, Gordan, and 
Alekseev and made use of their results, having filled them 
with quantum-mechanical content. 

As to isomer enumeration, the various attempts made 
in the 19th century to develop the Cayley approach any 
further failed. Only in 1930s the American scientists 
H. R. Henze and C. M. Blair derived the formulas for 
determination of the number of isomers in aliphatic com- 
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pounds of certain homological series. The solution of the 
problem discussed here was greatly promoted by the 
discovery of the so-called enumeration theorem of Polya 
in 1937*. 

Polya's theorem and benzene derivatives. Since the 
theorem is very complicated we cannot present it here, 
we shall analyse only one important partial case, the 
enumeration of isomers of substituted benzene deriva- 
tives. 

Table 6 shows the relationship between the symmetry 
group operations of topological graph of a benzene mole- 
cule and permutations of their vertices; the cyclic type of 
each permutation is also indicated. Using the data of 
this table, we can obtain the following expression simi- 
lar to (37): 

Z(C 6 H 6 ) = l/£ Y a{>a%, <», 

where g is the order of the group, i.e. the number of 
elements in the group. For the symmetry group of the 
topological graph of benzene g = 12. Therefore, in this 
case 

Z (C 3 H 6 ) = 1/12 {a\ + 2al + 2a\ + ?,a\a% + 4a*}. (39) 

To determine the number of isomers of substituted 
benzene derivatives (the substituting group is assumed 
to be monovalent) using Polya's theorem, a change of 
the variables must be performed in expression (39) 

i.e. 

a\-+{i + x)\ 
(l + *2)3, 

a^->-(l + a;) 2 (lra: 2 ) 2 - etc. 

After simple calculations we obtain the following 
enumeration polynomial 

1 + x + 3x 2 + Sx 3 + 3x 4 + x 5 + x*, 

* In a more partial formulation this theorem was given by 
mathematician D. Redfteld (1927) and independently by chemists 
A. Lann and J. Senior (1929). 
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Permutation of vertices of topological 
graph of benzene 
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in which the exponent of x is equal to the number of mono- 
valent groups bonded to the benzene ring, and the cor- 
responding coefficient indicates the number of isomers. 
Thus, the term 3x 2 indicates that there must be three 
disubstituted benzene derivatives (Fig. 51). 

Foresight of the Russian chemist. In connection with 
the problem considered above we would like to make one 
remark. In 1870-80s in chemical literature, primarily in 
Russian, the chemical structure theory was a controversial 
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subject provoking much discord. The leadings opponent 
of A. M. Butlerov was N. A. Menshutkin. D. I. Mende- 
leev also criticized the structure theory. They both pre- 
ferred the so-called substitution theory originating from 
the theory of types of Gerhardt. 

1 + x + Zx t + 3x 3 + Zx*+x s +x* 

oaaaaao 

» # « 
* * * 

0.00 

0.6.6 

Fig. 51 

According to Menshutkin, a molecule is an integral 
formation and its structure (in the sense adopted by 
Butlerov) cannot be determined by rigorous physical 
methods at the contemporary stage of scientific develop- 
ment, and therefore all hypotheses about the way the 
atoms are bonded should be discarded; the efforts should 
be concentrated on the study of similarity of properties 
of compounds. In Menshutkin's opinion, choosing from 
among the substances with similar properties the sim- 
plest ones, one may arrive atj the notion of chemical 
"type" of Gerhardt. "Classifying ajcompound with a type", 
pointed out Menshutkin, "implies that the same trans- 
formations can be produced with this substance as with 
the substance taken for the type." We shall not analyse 
here the strong and weak points of the) substitution the- 
ory; let us consider only the determination of the number 
of disubstituted benzene derivatives as proposed by 
Menshutkin. 

He proposed to consider, "without resorting to any 
structural] hypotheses", "the first product of benzene 
substitution in the order of both the direct and the re- 
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verse sequence", i.e. to enumerate the "fragments" 'of <i£H 
as follows: 

CH CH CH CH CH CX CX CH CH CH CH CH 
123456 654321 
direct sequence of carbon reverse sequence of carbon 
atom enumeration atom enumeration 



To determine "which fragments in the direct sequence 
correspond to the fragments in the reverse sequence", 
Menshutkin writes the foregoing formulas one under the 
other 

CH CH CH CH CH CX 
CX CH CH CH CH CH 

or, in the numerical notation: 



\ I 2 J 4 5\ 6 

6 15 4 3 2 /I 
i J 



As a result, the following permutation is obtained: 

/l 2 3 4 5\ 
P -{s 4 3 2 l) = (15 > (24) <3 > < 40 » 

Further arguments of Menshutkin can be easily inter- 
preted in terms of the permutation group theory. Here the 
number of cycles in the permutation (40) is equal to the 
number of disubstituted derivatives, i.e. 3 (ortho-, meta-, 
and para-), and the digits combined in one cycle indicate 
the numbers of equivalent positions in a monosubstituted 
derivative. 

Thus, although Menshutkin did not use explicitly the 
permutation groups (hardly being familiar with them at 
all), all of his arguments nevertheless are of "permuta- 
tion group" character. 

What was said in no way belittles the importance of 
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the Butlerov theory but only points out the profound 
rational content of the substitution theory proposed by 
Menshutkin* . 

Branches of a Topological Tree 

Up to now we dealt with counting isomers of various 
classes of chemical compounds. Not less important, how- 
ever, is to determine the topological indices characterizing 
each individual isomer. In the most profoundly studied 
case of structural isomers of alkanes the problem reduces 
to the determination of topological parameters defin- 
ing the extent of branching of their carbon skeleton. 
Unfortunately, many of the characteristics suggested for 
this purpose turned out to be insufficiently sensitive to 
replacing of one isomer by another. In 1977 D. Bonchev 
(Bulgaria) and N. Trinajstic (Yugoslavia) developed a meth- 
od of defining the extent of branching of a topological 
graph (more precisely, a tree) of alkanes; the method is 
based on the graph theory and the information theory. 

Information content. Originally, the term "information" 
denoted only the data exchanged by people and was of 
purely qualitative character. During the last 30 years the 
quantitative theory was elaborated, whose basic aspects 
will be briefly discussed below. 

Suppose a certain system comprises N elements. Then 
its information content is determined by Shannon's 
formula:** 

I = Nlog 2 N- |j JV,log 2 tf„ 

i=l 

where n is the number of different subsets of elements of 
the system combined according to some characteristic; N t 

* A more detailed discussion of this theory and the polemic 
between Butlerov and Menshutkin can be found in the article by 
I. S. Dmitriev and L. A. Kartseva "Chemical Structure Theory 
and Problem of Aromaticity" (Problems of history and methodology 
of chemistry. Issue 2. Ed. by R. B. Dobrotin, Leningrad University 
Press, 1978, in Russian). 

** The American scientist C. Shannon who proposed this for- 
mula in 1948 proceeded from the statistical interpretation of the 
concept of information. However, the elements can be selected not 
at random but according to some programme. The primary require- 
ment is the existence of a set of elements but not their random 
character. 
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is the number of elements in the i-th subset. The informa- 
tion content per one element of the system (/) is 

i 

where p t = NJN is the probability of finding an arbi- 
trarily chosen element of the system in the i-th subset. 

The information content is measured in special binary 
units, bits. According to W. R. Ashby information is 
a measure of variety of the system. Thus, the set consist- 
ing of three spheres of different colours possesses the 
variety of magnitude three. Every inhomogeneity of 
distribution of substance or energy in space is also expres- 
sed by the quantity / (or /). In the case of a completely 
homogeneous system (N { = N) I = 0, in the case of 
a completely inhomogeneous system (N t — 1) I = 
= I max = Nlog 2 N. 

In 1955 the US scientist N. Rashevsky introduced a con- 
cept of topological information of a graph (Itop)- His the- 
ory is based on subdivision of the set of N vertices of the 
graph into classes of topologically equivalent vertices. 
However, the quantity I top turned out to be unsuitable 
for characterization of branching of alkane's carbon skel- 
eton. 

Criterion of branching. The method of Bonchev and 
Trinajstic is based on the analysis of the matrix D 



Mi\0\l 2 3[ 

D,= ,2 1\0\1 21 I />,=> I 2 1 1 1 I D,= 



|3 2 1\0\1H 
ut_3__2__KO/ 




.p\l 2 3 4i\ /O 1 2 3 3\ 10 2 2 2 1 \ 




Fig. 52 



whose elements d kl represent the shortest distances be- 
tween the graph's vertices k and I (see Chapter 2 and 
Fig. 52). With increasing branching the distances in the 
graph become smaller. If the set containing iV 2 elements 
of the matrix D (N of which are zeros) is partitioned into 
groups as shown in Table 7, (m + 1) groups are obtained 
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where m = (d h i) max . The group containing all, zeros is 
not given in Table 7. The probability for a randomly 

Table 7 



Group No. 


1 


2 


3 


4 


Elements of 
matrix Di 


{4,4} 


{3,3,3,3} 


{2, 2,2, 2, 2, 2} 


{1,1,1,1,1,1,1,1} 



chosen matrix element d kl to be in the t-th group is 

Pi = 2k t /N\ 

where 2k t is the total number of elements in this group. 

Bonchev and Trinajstic introduced the concept of the 
information on distances within the graph (I D ) 

m 

I D = TV 2 log 2 N 2 — N log 2 N - 2 2fc/ log 2 2k, (41) 

i=l 

and 

m 

Ip=—jT log 2 x - 2 "FT lo S2 "77T • ( 42 ) 

i=l 

Noting that each element of the matrix D appears in 
its lower triangular submatrix the same number of times 
as in the upper submatrix (namely, k t times) (Fig. 52), 
Eqs. (41) and (42) can be replaced by the formulas; for 
one of the above-mentioned submatrices: 

,, N(N—l) , N(N—l) v 7 i 7 
Id= — ~2 — - iog2r4^ — ■— 2j*i. lo fif« & *». 

i 

*7' \i &i " i fc< 

1 D = — 2j N(N-i) l0gz N(N-j) 1 

1 2 , 2 

where — — is the total number of upper (or lower) 
off-diagonal elements in D. 

144 



For the topological tree of ra-pentane (N = 5) we obtain 
(Fig. 52): 

elements of matrix D x : 1, 2, 3, 4, 

values of k t : 4, 3, 2, 1, 
whence 

J ' D > = _ {lo 10+ lo iTJ + 

+ 15 log2 lo + lo log2 lo] = L84644 and = 18-4644. 

Table 8 



N 


Molecular graph 


li 




A 




R 754Q 


1.6 18 






6.0000 


1.732 


5 




to A R A A 

18 .4044 


1 71 9 
] . 1 O i 






15.2193 


1.848 






9.7095 


2.000 












32.2389 


1.802 


6 








28.6293 


1.902 






27.3842 


1.932 




o 1 1 o 


23.4840 


2.000 






22.5874 


2.074 







The quantities I'd and I' D characterize the inhomogene- 
ity of distribution of distances between vertices of a graph. 
Their values monotonically decrease as branching of 



10-01102 
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a molecule increases. Besides, the degree of blanching of 
a carbon skeleton can be compared with the greatest eigen- 
value of the topological matrix of a molecule (Table 8). 

The quantitative measures ^of molecular branching de- 
rived by Bonchev and Trinajstic are very important since 
they correlate with many physical and chemical proper- 
ties of compounds.* 

Transition State Topology 

Topological methods have found application not only 
in the investigations of the structure of a molecule but 
also in the study of mechanisms and rates of chemical 

Reaction: 

r .* - o - o 

dtene dienophile 

transition state: 




Fig. 53 



reactions; the topological theory of concerted reactions 
may serve as an example. 

Evans' rule. In 1938-39 the English chemist M. G. Evans, 
while studying the Diels-Alder reaction, assumed it 

* For more details see: Bonchev, D., N. Trinajstic, /. Chem. 
Phys. 67, 10, 4517 (1977). 
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to be a synchronous process in the course of which six 2p 
carbon AOs (four of diene and two of dienophile) overlap 
to form a cyclic transition state which is topologically 
equivalent to benzene * (Fig. 53) and therefore stabilized. 
In its turn, the process of ethylene dimerization to form 
cyclobutane must proceed through the anti-aromatic 
(destabilized) transition state topologically equivalent 

Reaction : 

M - n - □ 



Transition state: 




Fig. 54 



to cyclobutadiene (Fig. 54); it follows that this reaction 
is thermally forbidden. 

Evans' ideas were developed in 1952 by M. Dewar and in 
1966 by H. E. Zimmerman. The synchronous reactions 
whose transition state can be represented as a cycle to 
which each constituent atom donates one AO are referred 
to as pericyclic. 

Anti-Hiickel systems. In 1964 E. Heilbronner considered 
an unusual type of cyclic at-electron systems with polyene 
ribbon twisted an odd number of times. The topology of 



* The imaginary connection of overlapping lobes of 2p-AO 
transition state results in a curve which is topologically equivalent 
(homeomorphic) to the curve obtained in the similar analysis of 
a benzene molecule. The difference between a and jt overlaps of AOs 
is not taken into consideration. 



10* 
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such a molecule resembles the Mobius band, (Fig. 55), 
a surface first described by the German mathematician 
J. B. Listing in 1861 and later, in 1863-65, investigated 
by A. Mobius in full detail. The Mobius band cannot be 
painted in two colours, one inside and another outside, 
because it is a one-sided surface. The cyclic n-electron 
systems possessing the Mobius band topology are referred 
to as anti-Hiickel, or Mobius systems. In such systems there 
are an odd number of spots where the AOs of the neigh- 
bouring! atoms overlap "out nf phase". According to 

AO phase shift 

g. 55 

Dewar, "although none of such systems has been obtained 
so far, it is of great significance to develop a conception 
to substantiate it. The system in which the AO phase 
shift is inevitable differs topologically from the system 
in which such a shift can be avoided... . This difference 
depends only on the overlap topology of AOs involved in 
the formation of a delocalized system and is absolutely 
independent of the MOs that can be composed from 
those AOs". 

Anti-Hiickel aromaticity. The aromaticity rules for 
anti-Hiickel systems are opposite to those for conventional 
(Huckel) systems (Table 9). 



Table 9 



Number of 
n-electrons 


Topological type 
of system 


Character of system 


4JV + 2 


Huckel 


Aromatic 




Anti-Hiickel 


Anti-aromatic 


AN 


Huckel 


ditto 




Anti-Hiickel 


Aromatic 
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Zimmerman's diagrams. For anti-Hiickel systems thee 
MO energies are defined by the following formula: 

OQ [~2n (fc + 1) | 
e ft = a-2pcos [ — 7^ j 

(fc = 0, 1, 2, ...,7V-1). 

Accordingly, the mnemonic Frost rule is replaced by 
the Zimmerman rule for these systems. The latter differs 
from the former by the position of the inscribed regular 
polygon (Fig. 56, cf. Fig. 5). 



— +2 




+0.618 

0.000 

1.618 

+1.732 




Evans- Dewar-Zimmerman criterion. In accordance with 
this criterion the pericyclic reactions proceed through 
the cyclic transition state of the Hiickel topological type 
if it contains (47V + 2) n -electrons, and through the state 
with anti-Hiickel (Mobius) topology if it contains 47V 
ji -electrons. 

To determine the topological type of a transition state, 
one should "walk" around the cycle and count how many 
times the AO changes its sign*. This sign changes occur- 



* Surely, the cycle can contain not only p AOs. 
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ring while we move along the same AO are ignored. The 
Hiickel system has an even number of sign changes 
(including zero), the anti-Hiickel one has an odd number. 

As an example, let us consider the electrocyclic reac 
tion of cyclobutene to butadiene transformation. The 
cyclobutene ring can break up either in a conrotatory 




Fig. 57 



way (Fig. 57a) or in a disrotatory way (Fig. 576). In 
both con- and disrotatory transition states the hybrid AOs 
of carbon methylene atoms interact not only between 
themselves but also with p„ AO of the C = C bond, i.e. 
each carbon atom of the system possesses an AO interact- 
ing with AOs of the neighbouring atoms in the ring. Here 
the four-electron disrotatory transition state relates to 
the Hiickel type and is therefore anti-aromatic whereas 
the conrotatory state associated with the anti-Hiickel 
topology is aromatic; consequently, the indicated reac- 
tion proceeds in a conrotatory way. 

Symmetry or topology? For the first time synchronous 
reactions were theoretically studied by R. B. Woodward 
and R. Hoffmann on the basis of the orbital symmetry 
conservation principle formulated by them. This princi- 
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pie procAumed initially in the fashion of the Sermon on the 
Mount wad subsequently subjected to a detailed quantum- 
mechanical analysis. H. Metiu, J. Ross, and G. M. White- 
sides* demonstrated that the Woodward-Hoffmann rules 
emerged from the most simplified treatment in which: 

(1) the spin-orbital coupling, the virtual transitions 
between excited electronic states of reactants and products 
and the non-Born-Oppenheimer terms are ignored; 

(2) the Franck-Condon approximation is used, and 

(3) the electronic wave functions are determined by the 
HMO method. Therefore, a more complete calculation 
may, in principle, lead to predictions different from those 
of the Woodward-Hoffmann rules. 

In addition, Dewar demonstrated that the classification 
of pericyclic reactions has nothing to do with symmetry. 
The nature of the reaction is defined by the AO overlap 
topology in a pericyclic transition state and not by tbe 
MO symmetry. 



We have described various applications of the graph 
theory in chemistry. Surely, much was left outside this 
book, in particular such problems as topology of molecu- 
lar potential energy surface, topological aspects of ther- 
modynamics of heterogeneous equilibrium systems, the 
graph theory and coding of chemical information. How- 
ever, even what was said suggests the fundamental impor- 
tance of topological ideas and methods. The primary areas 
of chemical science in which the graph-theoretical and 
topological methods are employed are shown in Fig. 58. 

The topological ideas play an important part in ste- 
reochemistry. Here is one fairly indicative example. In 
the latest chemical literature the problems associated 
with the concepts of stereochemical configuration, con- 
formation, and chirality have been extensively discussed. 
In particular, the'concept of configuration was shown to 
be'applicable'only to chiral molecules, i.e. to such mole- 



* H. Metiu, J. Ross, G. M. Whitesides. A Basis for Orbital 
Symmetry Rules.— Angew. Chem. Intern. Edition in English, 18, 5, 
377 (1979). 



Conclusion 




cules that cannot be superimposed upon their mirror im- 
ages by translation and rotation. Hence, various types of 
molecular deformations proceeding both with a change 
in configuration and with the preservation of the latter 
were considered. Of special interest is the first type of 




Fig. 58 



transformation going either via chiral (e. g. Berry's pseu- 
dorotation for a trigonal bipyramid, see Fig. 47) or achiral 
(e. g. in the diagonal twist) intermediate. In either case 
the topological pattern) of a molecule changes. When 
considering such intramolecular motions it should be 
taken into account that a molecule may possess some 
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singulai points; this fact makes it necessary to introduce 
certain c6rrections into the strictly geometrical consid- 
eration. For example, let us look at a molecule of symme- 
try T d . As distinct from the tetrahedron treated in a strict 
geometrical sense (and we remind the reader that in that 
sense all tetrahedrons are homeomorphic), a tetrahedral 
molecule, as a rule, possesses a singular point, the point 
at which the central atom is located. When such a tetra- 
hedral structure is transformed into a disphenoid one 
and then into a tetragonal pyramid, this singular point 
(the central atom) first comes out to the surface of the 
disphenoid and then becomes the pyramid's vertex. This 
and other similar examples have brought about the con- 
cept of molecular topological form which is understood as 
a topological pattern in combination with singular points 
of the molecule. The molecular topological form was 
shown to be unique for every molecule*. If the molecular 
topological form is retained in the process of a conforma- 
tional transition, the conformations associated with that 
transition are of the same configuration type. In other 
words, the configuration remains unchanged for a certain 
set of molecular conformations. If we now take into 
account that the concept of configuration (the absolute 
configuration is meant) is meaningful only for chiral 
molecules, the following definition can be given: the 
stereochemical (absolute) configuration is an inherent 
invariant of a chiral molecular topological form. 

Such an approach turned out to be very fruitful in the- 
oretical stereochemistry. In particular, it helps to estab- 
lish a topological resemblance between various chemical 
structures. For example, it can be shown that the confor- 
mational transition between the eclipsed and staggered 
conformations of an ethane molecule is similar to the 
trigonal twist in an octahedron resulting in a trigonal 
prism formation. In these extreme conformations the 
molecular topological form of ethane reduces to a trigo- 
nal prism and a trigonal antiprism possessing two internal 
singular points occupied by carbon atoms. 



* The readers interested in problems of chemical topology and 
application of topological ideas to stereochemistry are advised to 
look through the remarkable book by V. I. Sokolov Introduction 
to Theoretical Stereochemistry, Moscow, 1979 (in Russian). 
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At present, in stereochemistry, as in other, areas of 
chemistry, ever increasingly important for traditional 
chemical problems become such unusual (but on the 
second thought, quite natural) mathematical disciplines 
as the set theory, topology, the group theory, mathemati- 
cal logic, the information theory, etc. Here we witness 
the nonquantitative mathematization of knowledge, so 
characteristic for modern natural sciences. 
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